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Summary. — Several topics on the implementation of spin qubits in quantum dots 
are reviewed. We first provide an introduction to the standard model of quantum 
computing and the basic criteria for its realization. Other alternative formulations 
such as measurement-based and adiabatic quantum computing are briefly discussed. 
We then focus on spin qubits in single and double GaAs electron quantum dots 
and review recent experimental achievements with respect to initialization, coherent 
manipulation and readout of the spin states. We extensively discuss the problem 
of decoherence in this system, with particular emphasis on its theoretical treatment 
and possible ways to overcome it. 



1. — Introduction 



It was in the 1980s, when the idea of exploiting quantum degrees of freedom for 
information processing was envisioned. The central question at the time was whether and 
how it was possible to simulate (efhciently) any finite physical system with a man-made 
machine. Deutsch [1] argued that such a simulation is not possible perfectly within the 
classical computational framework that had been developed for decades. He suggested. 
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together with other researchers such as Fcynman [2, 3], that the universal computing 
machine should be of quantum nature, i.e., a quantum computer. 

Around the same time, developments in two different areas of research and industry 
took a tremendous influence on the advent of quantum computing. On the one hand, 
it was experimentally confirmed [4] that Nature indeed does possess some peculiar non- 
local aspects which were heavily debated since the early days of quantum mechanics [5] . 
Schrodinger [6] coined the term 'entanglement', comprising the apparent possibility for 
faraway parties to observe highly correlated measurement results as a consequence of 
the global and instantaneous collapse of the wave function according to the Copenhagen 
interpretation of quantum mechanics. The existence of entanglement is crucial for many 
quantum computations. On the other hand, the booming computer industry led to 
major progress in semiconductor and laser technology, a prerequisite for the possibility 
to fabricate, address and manipulate single quantum systems, as needed in a quantum 
computer. 

As the emerging fields of quantum information and nanotechnology inspired and mo- 
tivated each other in various ways, and are still doing so today more than ever, many 
interesting results have been obtained so far, some of them we arc; about to review in this 
work. While the theories of quantum complexity and cntanglcnicint are being cstablish(xl 
(a process which is far from being complete) and fast quantum algorithms for classically 
difficult problems have been discovered, the control and manipulation of single quantum 
systems is now experimental reality. There are various systems that may be employed 
as qubits in a quantum computer, i.e., the basic unit of quantum information. Here, we 
will focus on the idea of using the spins of electrons confined in quantum dots, a proposal 
made in 1997 [7]. Most knowledge for realizing qubits in the fashion of the spin-qubit 
proposal of rcf. [7] has so far been obtained for quantum dots formed in a two-dimensional 
electron gas at the interface of a GaAs/AlGaAs heterostructure. Therefore, this is the 
main type of quantum dot we will turn our attention to. We would nevertheless like to 
mention that many other systems are investigated intensively at present, such as car- 
bon nanotubes [8], nanowires [9], molecular magnets [10, 11, 12, 13], quantum dots in 
graphene [14], and nitrogen- vacancy centers in diamond [15, 16, 17, 18]. 

The idea behind this work is to review quantum computing starting at its very roots, 
and ending at the current state of knowledge of one of its many branches, here being the 
spin-qubit proposal of ref. [7] for electron spins in GaAs quantum dots. We approach 
the subject in the second section by beginning with classical computing and complexity 
theory, mainly to sec the concepts that inspired the so-called standard model of quantum 
computing discussed later, and to encounter examples of problems that are assumed to 
be too involved to ever be solved in reasonable time on a classical computer. We then 
review a set of requirements that every quantum computing proposal should fulfill in 
order to be usable on a large scale [19]. After discussing the actual spin-qubit proposal 
of ref. [7] , we end the chapter with some recent alternative models of quantum computing 
and a section about entanglement measures and their numerical evaluation. 

While the second section is quite general, we focus in the third section on spin ma- 
nipulation in GaAs quantum dots. We discuss the latest experimental achievements and 
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we will see that it recently became possible to define, initialize, manipulate, and read out 
single electron spin qubits with an already quite remarkable rate of success. 

Although anticipated to some extent in the third section, the last section intensively 
discusscis the various mechanisms of dccohcrcncc and possible ways to reduce its effects, 
which is a necessary prerequisite for large-scale quantum computing. In greater detail, 
we investigate the role of the spin-orbit and the hyperfine interaction. We show how these 
mechanisms cause relaxation of the electron spin, but also how they can be exploited to 
perform such beneficial tasks as all-electrical single spin manipulation or full polarization 
of the quantum dot's nuclear bath. 



2. — Quantum computing in a nutshell 

2'1. Classical computers and complexity theory. - Computers are devices that solve 
problems in an algorithmic fashion. The Church- Turing hypothesis [20, 21] claims that 
every function which we would naturally regard as being computable by an algorithm 
{i.e., a procedure that solves a given problem in finite time) can be computed by the 
universal Turing machine (i.e., a hypothetical, mathematically formalized, programmable 
discrete machine). Fortunately, the rather cumbersome notion of a Turing machine 
turns out to be computationally equivalent to a more human-friendly description of 
algorithms, namely the circuit representation. The basic unit of information, the bit, is a 
physical system that can be in exactly one out of two states, usually denoted by and 1. 
Information (or data) is stored in binary form using many bits which are represented as 
lines in the circuit. The data is processed by consecutively applying logical gates to the 
bits, depicted symbolically as elements acting between the lines in the circuit, such that 
the desired algorithm is performed and the solution to the problem is stored in the bits. 
Examples of such gates include the 1-bit not gate which flips the state of a bit from to 
1 and vice versa, and the 2-bit XOR gate that outputs 1 if and only if exactly one of the 
input bits is 1 (see fig. 1). An important result here is the fact that there exists a set of 
gates with which one can implement any algorithm in a circuit, provided one can freely 
distribute and copy information. For classical information, the latter constraint is trivial 
and the so-called universal set of gates consists in fact of only the nand gate (yielding 
if and only if both input bits are 1). 

Having set the playground to implement algorithms tackling computational problems, 
one of the most profound questions one can ask is the following: What is the most 
efficient algorithm to solve a particular task? It is the field of computational complexity 
theory [22] that deals with such kind of issues. We will briefly discuss the most important 
complexity classes, focussing on time rather than space complexity. Imagining that every 
gate in a circuit requires a finite execution time, one can study the total time required to 
run an algorithm as a function of the input or problem size, e.g., the number of input bits 
n. A procedure is called efficient (or tractable) if its running time is upper-bounded by 
some polynomial in n. Roughly speaking, the complexity class P consists of all problems 
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Fig. 1. - a) Classical XOR gate, b) Its quantum analog, i.e., the CNOT gate. The truth table of 
the XOR gate is shown in the lower part of a) and simply comes from adding the two input bits 
modulo 2 (this is the '©' operation). By replacing with |0) and 1 with |1) one obtains the 
truth table for the target qubit (lower wire) in CNOT if the computational basis states are used 
as input. However, the CNOT gate is also the source of entanglement (see also subsect. 2'7) in 
the standard model of quantum computing, as indicated in the lower part of b) by the fact that, 
generally, two states cannot be written as a product state ® jV'2) anymore after CNOT has 
been applied. 



known to have efficient algorithms ("'^). Sorting or searching lists are examples thereof. 

However, there is also a vast amount of problems solved by algorithms that are not 
efficient (typically exponential in n), but once a possible solution has been proposed, it 
can be efficiently checked for its validity. The class of these types of problems is called 
NP. An important example thereof is integer factorization: Given an n-bit integer, the 
best known classical algorithm to find its prime factors is exponential in n, but given a 
proposed factorization one can quickly check whether it is correct by doing the required 
multiplications. The observation that there are problems that can be solved efficiently, 
and others, for which the best known solution is still worse than polynomial, cumulates 
in the famous P = NP? question: It is clear that P is a subset of NP, but is this inclusion 
strict? In other words, have we just not yet discovered efficient algorithms for supposedly 
'hard' problems, or is there something fundamental within these problems that prevents 
us from finding such? This puzzling question has not been answered yet with a formal 
proof, but it is widely assumed that P and NP are not equal. 

Another important complexity class is NP-complete, a subset of NP consisting of all 
problems that are at least as hard as all other problems in NP. This means that every 
problem in NP can be cast into an instance of a problem in NP-complete in polynomial 
time. With this transformation being efficient, a polynomial-time solution to any one of 



(^) We omit here the strict definitions of complexity classes in terms of decision problems and 
formal languages. The interested reader is referred to, e.g., ref. [22]. 



Quantum Computing with Electron Spins in Quantum Dots 



5 



the problems in NP-complete would render the whole class NP tractable and it would 
follow that P = NP. Integer factorization is widely suspected, but not proven, to be both 
outside of P and NP-complete. 

Finally, we remark that there is a complexity class of high practical interest, which 
we include here merely because it has a quantum analog that we will encounter in sub- 
sect. 2'3. Decision problems in the class BPP (for bounded-error probabilistic time) have 
efficient algorithms that arc allowed to make random choices ('coin flipping') during the 
computation and yield the correct answer with probability p > 1/2, and a wrong solution 
with probability 1 — p. The choice of p is essentially arbitrary, since the Chernoff bound 
[23] guarantees that the error probability in a majority vote drops exponentially with the 
number of repetitive executions of the algorithm. Typically, however, one finds p = 2/3 
or p = 3/4 in the literature. 

2"2. The standard m,odel of quantum computing . - In a seminal work by Dcutsch 
[1], he proposes to strengthen the Church- Turing hypothesis into a "manifestly physical 
and unambiguous" form. His Church- Turing principle reads: "Every finitely realizable 
physical system can be perfectly simulated by a universal model computing machine 
operating by finite means", arguing that "it would surely be hard to regard a function 
'naturally' as computable if it could not be computed in Nature, and conversely". It 
is further shown that the universal Turing machine does not fulfill this principle, while 
the 'universal quantum computer', proposed in the same work, is compatible with the 
principle. A maybe less theoretic reasoning for studying computing machines operat- 
ing in the quantum regime is the mere fact that classical computers are governed by 
Newtonian mechanics, being valid only in a limiting case of the underlying quantum 
theory [19]. Quantum computers must therefore have at least the same, if not a greater, 
computational power than classical computers. 

Historically understandable, the now so-called standard model of quantum computing 
[19; 24] follows closely the circuit model discussed earlier (for alternative proposals of 
quantum computing, see subsect. 2'6). The basic unit of information is the quhit, being 
a two-level quantum system with basis states usually denoted by |0) and |1) according to 
its classical counterpart. Qubits are displayed as lines in the circuit to which quantum 
gates are applied successively, thereby performing the computation. The final result is 
obtained as a readout (a measurement) of the qubits in the computer's final state. 

The following features distinguish quantum from classical computing. Firstly, quan- 
tum states cannot be copied perfectly (this is the no-cloning theorem [25]). Secondly, 
the nature of a quantum computer is ultimately an analog one. Quantum gates operate 
on amplitudes such as a and /? in the state = a\(pi) + (3\(f)2). Apart from being prop- 
erly normalized, ampHtudes are arbitrary complex numbers and as such analog. Thirdly, 
qubits can be in superpositions and may form intricate entangled states. This fact is 
heavily exploited in existing quantum algorithms and is essentially the key ingredient to 
the majority of them. These first three points all contribute to the drawback that error 
correction, a necessity in the presence of imperfect gates and decohcrence, is a non-trivial 
thing to do. Quantum error correction, however, turns out to be possible if the error 
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probability per gate is smaller than some finite value for all gates (sec criterion 3 in 
subsect. 2'4). And lastly, quantum gates must be time-reversal, i.e., unitary operators, 
in accordance with basic principles of quantum mechanics. Apart from the latter, there 
are no other constraints imposed on quantum gates. 

Analogous to the classical case, and quite remarkably, there exist finite sets of gates 
which can be used to approximate any unitary evolution {i.e., the computation) of the 
quantum machine to arbitrary precision [26]. An example thereof is the universal set 

[24] consisting of the two single-qubit gates = ^ j and T = 

and the two-qubit CNOT gate which performs the operation \x)\y) \x)\{x + y) mod 2), 
x,y € {0, 1}. The gates are represented as usual in the standard computational basis, i.e., 
|0) = (1,0)-^ and |1) — (0, 1)-^. The H- and T- gates are used to approximate arbitrary 
single-qubit rotations e*"e~'^"'°'/^, where n is a real unit vector and tr is the vector of 
Pauli matrices. If arbitrary rotations arc available on their own, one can together with 
the CNOT gate implement any unitary evolution exactly. 

2'3. Quantum algorithms and quantum complexity. - Having the new quantum ma- 
chine at hand, what is its computational power, and how does it perform compared with 
classical computers? Such questions have given birth to the new field of quantum com- 
plexity theory [27, 28] , resulting in a plethora of new quantum complexity classes along 
with the goal of understanding their relations both between each other, and to classical 
complexity classes. This is an active field of research and the goal mentioned before is far 
from being reached. Instead of going into the details of the theory, which would be out 
of the scope of this work, we sketch an overview of the subject using the most prominent 
examples. 

One of the first algorithms demonstrating that quantum computers may be able to 
drastically outperform classical computers is the Deutsch-Jozsa algorithm [29, 30]. Let 
/ be a function from n bits to one bit that can either be constant or balanced, the 
latter meaning that there are exactly 2"~^ unknown input strings yielding the output 0. 
The function is supposed to be implemented in an oracle, i.e., a black box without any 
further internal specification. The Deutsch-Jozsa algorithm can determine the function's 
type {i.e., constant or balanced) with probability 1 by querying the oracle just once. 
A deterministic classical algorithm requires 2"~^ -|- 1 queries in the worst case, as one 
may coincidentally pick all the 2"~^ input strings that yield the same output (^). The 
Deutsch-Jozsa algorithm obtains its power from bringing a register of n qubits into an 
equally weighted superposition of all possible bit strings using a Hadamard transform 
which is a special case of a so-called quantum Fourier transform [24]. Using this state as 
the input for the oracle, the function can be evaluated simultaneously on all strings in the 
superposition with just one query. This functionality in quantum computing is referred to 




(^) In practice, however, one would use a randomized algorithm requiring a constant number of 
k queries and returning the wrong result with arbitrary low probability (dropping exponentially 
with k). 
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as quantum parallelism. Without going into further details, the Deutsch- Jozsa algorithm 
manages to output the correct result with certainty using these techniques. However, one 
has to be cautious not to get the impression that one can calculate and obtain all function 
values of an arbitrary function with just one query of an oracle. In general, one is left with 
a superposition of results collapsing upon measurement and yielding just one function 
value. The Deutsch-Jozsa algorithm is cleverly designed to work with certainty for a 
particular type of problem. These ideas arc not straightforwardly adapted to problems 
involving other kinds of functions. The constant/balanced problem is an example from 
the class EQP (exact quantum polynomial-time) [27], denoting all problems solved in 
polynomial time by a quantum algorithm with success probability equal to 1. EQP is the 
analog to the classical complexity class P. 

On the other hand, the set of problems having tractable algorithms on a quantum 
computer with error probability smaller than p = 1/3 is denoted by BQP (bounded- 
error quantum polynomial-time) [27], its classical counterpart being BPP. It contains 
more interesting problems whose solutions are of greater relevance than the Deutsch- 
Jozsa algorithm. Famous examples thereof are Grover's [31] and Shor's [32] algorithms. 
Grover's algorithm searches an unsorted database in time upper-bound by a function 
proportional to \/]V and has a probability of failure scaling as where N is the 

number of database entries. Although the speed-up with respect to the classical database 
search which runs in time proportional to N is not exponential, it is still of great benefit 
especially for large N . Grover's algorithm could be used to speed up brute-force search 
attempts for finding solutions to computationally hard problems. However, probably the 
most prominent problem in BQP is integer factorization. Shor [32] has shown that a 
quantum computer can factor an integer in polynomial time, a task which is supposed to 
be exponentially hard on a classical computer. These kind of discoveries have resulted in 
a tremendous increase of interest in quantum computing. For example, Shor's algorithm 
can be used to break present-day public-key cryptosystems who rely on the hardness of 
factorizing integers being a product of two large prime numbers (such as RSA [33]). 

A quantum complexity class that has recently called attention is QMA [34, 35] (quan- 
tum Merlin- Arthur), which can be seen as a quantum analog of NP. Equivalently to the 
definition in subsect. 2"1, NP can be characterized in terms of decision problems (hav- 
ing either 'y^s' or 'no' as answer): NP contains all problems for which 'yes'-instances 
are supplied with a proof that can be checked in polynomial time by a deterministic 
verifier(^). For quantum computers, deterministic verifiers are however not meaningful. 
QMA is thus defined probabilistically: A decision problem belongs to QMA, if for every 
instance x there exists an efficient {i.e., polynomial in x) description of a quantum cir- 
cuit (the verifier), such that for every 'yes'-instance x there exists a proof lip) with 



(^) For example, a decision version of integer factorization would be stated as: Given integers 
A'' and a (user-specified) x < N, is there an integer 1 < d < a; such that d divides N? The 
answer 'yes' provided with a suitable d can be verified in polynomial time, simply by checking 
whether 1 < d < x and d divides N. One then says that x provides a 'yes'-instance to the 
problem with proof d. 
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p{Qx accepts IV')) > 2/3, and for every 'no'-instance x p{Qx accepts |x)) < 1/3 holds for 
all input states |x)- Here, p{Qx accepts \^p)) denotes the probability to measure, e.g., |0) 
(if this is defined to mean 'accept') as the output of a quantum computation described 
by the circuit and started with the initial state \ip). 

An important problem known to be in QMA is fc-LOCAL Hamiltonian which is 
specified by the following decision problem: Given a Hamiltonian H acting on n qubits 
with interac;tions that do not involve more than k particles (/c-body interactions, where k 
is a constant independent of n) and two real numbers a and b, such that 6— a > 1 /poly (n) . 
Is the ground state energy of H smaller than a ('yes'), or are all energies larger than b 
('no')(^)? A 'yes' instance can be verified by providing an eigenstate with energy smaller 
than a. Furthermore, one can show that polynomial verifiers can be constructed that 
accept (reject) 'yes' instances ('no' instances) with sufficient probability. It is now known 
that /c-LOCAL Hamiltonian is QMA-complete for > 2 [36], meaning the following: 
Given an instance x of any problem Q in QMA, one can find (in time poly(n)) an instance 
of fc-LOCAL Hamiltonian (by constructing a fc-local H and specifying properly chosen 
parameters a and b), such that, if x is a 'yes' instance of Q, the ground state energy of H 
is smaller than a, and if a; is a 'no' instance, the smallest eigenvalue of H is larger than 
b. fc-LOCAL Hamiltonian (for fc > 2) is as hard as any other problem in QMA. This 
'hardness' suggests that calculating the ground state energy, and possibly other ground 
state properties, is intractable even on a quantum computer. The study of fc-LOCAL 
Hamiltonian is also important in the context of adiabatic quantum computing, see 
subsubsect. 2 '6. 2. 

2"4. General criteria for scalable quantum computing. - In this section, we review 
the five DiVinccnzo criteria [19] for the physical implementation of quantum computing. 
These are the most fundamental requirements any proposal for a quantum computer 
must fulfill in order to work with an arbitrary number of qubits. Starting from sect. 3, 
we will examine the experimental and theoretical progress toward realizing these criteria 
for the spin-qubit proposal of ref. [7] (see the next section). 

1. A scalable physical system with well characterized qubits. 
We have already defined the notion of a qubit as simply being a two-level quantum sys- 
tem. In this review, we will focus solely on the electron spin in quantum dots. The 
word 'scalable' plays an important role: Even if current fundamental experiments are 
performed with only few qubits, they must at least in principle be preparable or man- 
ufacturable in large numbers, since only in this case interesting and useful quantum 
computations can be performed. 

A qubit must also be 'well characterized' in the sense that one has a good theoretical 
description not only of the qubit itself (in terms of an internal Hamiltonian, accurate 
knowledge of all physical parameters, etc.), but also of all relevant mechanisms that 



(*) Note that this is a promised problem: It is guaranteed that either of the two cases will 
occur, and we are not interested in Hamiltonians that have energies between a and b. 
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couple qubits among each other and to the environment. On the one hand, this is 
necessary to explore the possibilities of manipulating qubits and letting them interact, 
but also, on the other hand, in order to understand and fight the various forms of 
decoherence a qubit may suffer from (sec sect. 4). 

2. The ability to initialize the state of the qubits to a simple fiducial state. 

It is clear that every computation needs to be started in an initially known state such as 
|000 . . .). But this is not the end of the story. Having a fast initialization mechanism at 
hand is crucial for quantum error correction (see next criterion), typically requiring large 
amounts of ancillary rjubits in known initial states in order to perform its job properly. 
If a fast zeroing of qubits is not possible, i.e., if the initialization time is long compared 
to gate operation times, then ref. [19] proposes to equip the quantum computer with 
"some kind of 'qubit conveyor belt', on which qubits in need of initialization are carried 
away from the region in which active computation is taking place, initialized while on 
the 'belt', then brought back to the active place after the initialization is finished." 

In spin qubits, initialization could be achieved by either forcing the spins to align 
with a strong externally applied magnetic field, or by performing a measurement on 
the dot followed by a subsequent rotation of the state depending on the measurement 
outcome. The first approach is somewhat problematic, since natural thermalization times 
are always longer than the decoherence time which itself needs to be much longer than 
gate operation times (see the next criterion). In this case, a 'conveyor belt' scheme would 
be required. The second possibility of measurement and rotation depends on the specific 
setup examined, but initialization times might in principle be much shorter than natural 
relaxation times. See subsects. 2'5 and 3"2 for more information and recent experimental 
achievements. 

3. Long relevant decoherence times, much longer than the gate operation time. 

Due to the coupling of qubits to their environment in a thermodynamically irreversible 
way, quantum coherence is lost. In other words, quantum states in contact with the 
outside world ultimately evolve into fully mixed states. Decoherence is the answer to 
why the macroscopic world looks classical. In GaAs quantum dots where electron spins 
are used as qubits, the most important mechanisms of decoherence are the spin-orbit and 
the hyperfine interaction, see subsect. 3'3 and sect. 4. 

If it were not for quantum error correction, the duration of a quantum computation 
would eventually be determined by the shortest decoherence time in the setup. This 
would render longer and more complex computations impossible. By encoding informa- 
tion not directly into single qubits, but rather into 'logical qubits' consisting of several 
single qubits, a certain amount of errors due to decoherence and imperfect gates may 
be corrected, depending on what kind of code is used. There is however still a limit 
on how faulty elementary gates are allowed to be: The accuracy threshold theorem [37] 
states that error correction is possible if the error probability per gate is smaller than a 
certain threshold. This threshold comes about the fact that encoding, verification, and 
correction steps require an additional overhead of quantum gates, which introduces new 
possible sources of errors. Error correction is thus only meaningful if encoding reduces 
the error probability of an encoded operation compared with its original ('unencoded') 
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counterpart, despite the fact that more elementary gates are required. By employing 
the technique of code concatenation [37], i.e., encoding logical qubits recursively, where 
using different codes per concatenation level is allowed, the effective error on the top 
level of concatenation can be made arbitrarily small. The threshold value depends on 
the error models studied and on the details of the codes considered. Typical values are in 
the range of 10~^ to 10~^ [37, 38], implying that decoherence times must be a thousand 
to a hundred thousand times longer than gate operation times. 

4- A "universal" set of quantum gates. 
We have mentioned in subsect. 2'2 that generic quantum computing is possible in the 
standard model if certain one- and two-qubit gates are available. The single qubit gates 
may be either implemented directly, or can be approximated to arbitrary precision using 
a finite set of gates. The only necessary two-qubit gate is the controUed-NOT gate [26]. 
If these gates are for some reason not implementable directly {i.e., in the sense that 
there are no Hamiltonians that can be switched on which perform exactly the desired 
gate operations), then a set capable of synthesizing them needs to be present. This 
is, e.g., the case in the spin-qubit proposal of ref. [7], where the only controllable two- 
qubit interaction is the exchange coupling between neighboring spins. The CNOT gate 
can however be implemented using a series of one-qubit operations and the exchange 
interaction alone. The details and requirements for this to work are discussed in the next 
section. 

It is worth pointing out that one also needs to be able to execute quantum gates 
in parallel in order for error correction to work. This does however not pose a major 
drawback for solid-state systems [39] , where usually only two-body nearest-neighbor in- 
teractions are realizable. It is also important to note the fact that faulty gates may 
introduce systematic or random errors in a calculation. This can be viewed as a source 
of decoherence and can therefore be overcome by means of quantum error correction if 
the error rate is sufficiently small. The same threshold values as discussed in the previous 
criterion hold in this case. 

5. A qubit- specific measurement capability. 
Measuring qubits without disturbing the rest of the quantum computer is required in 
the verification steps of quantum error correction and, not remarkably, in order to reveal 
the outcome of a computation. If the measurement procedure does not discard qubits 
(which could be the case, e.g., for spin-dependent tunneling of electrons out of a quantum 
dot) it may be used in the initialization step (see 2nd criterion). If it is, additionally, 
fast enough, it may also be useful for quantum error correction. A measurement is said 
to have 100% quantum cfRciency if it yields, performed on a state p = p\0){0\ + (1 — 
p)|l)(l| -ha|0)(l| -ha*|l)(0|, the outcome "0" with probability p and "1" with probability 
(1 — p) independent of a, the states of neighboring qubits, or any other parameters of 
the system. Real measurements cannot have perfect quantum efficiency. But this is also 
not required since one can, e.g., rerun the computation several times. 

2'5. The Loss-DiVincenzo proposal. In this section, we review the spin-qubit pro- 
posal of ref. [7] for universal scalable quantum computing. Here, the physical system 




high-g layer quantum well 

Fig. 2. - An array of quantum dot qubits realized by laterally confining electrons in a two dimen- 
sional electron gas formed at the interface of a heterostructure. The confinement is achieved 
electrostatically by applying voltages to the metallic top gates. Interaction is generally sup- 
pressed (as for the two qubits on the left) but may be turned on to realize two-qubit operations 
by lowering inter-dot gates (as for the two qubits on the right). Single spin rotations may be 
achieved by dragging electrons down (by changing back gate voltages) to a region where the 
Zeeman splitting in the presence of the external static magnetic field B± changes due to magne- 
tization or an inhomogeneous g-factor present in that layer. A resonant magnetic ac pulse _B||'^ 
can then be used to rotate the spin under consideration, while leaving all other qubits unaffected 
due to the off-resonant Zeeman splitting (ESR). All-electrical single spin manipulation may be 
realized in the presence of spin-orbit interaction by applying ac electric pulses directly via the 
gates (EDSR). See subsect. 3'4 for more details. 



representing a qubit is given by the localized spin state of one electron, and the com- 
putational basis states |0) and |1) are identified with the two spin states |t) and ||), 
respectively. In general, the considerations discussed in ref. [7] are applicable to elec- 
trons confined to any structure, such as, e.g., atoms, defects, or molecules. However, 
the original proposal focuses on electrons localized in electrically gated semiconductor 
quantum dots. The relevance of such systems has become clearer in recent years, where 
remarkable progress in the fabrication and control of single and double GaAs quantum 
dots has been made (see, e.g., ref. [40] for a recent experimental review). We postpone 
the discussion of experimental achievements with respect to satisfying the DiVincenzo 
criteria to sect. 3. 

Scalability in the proposal of ref. [7] is due to the availability of local gating. Gat- 
ing operations are realized through the exchange coupling (see below), which can be 
tuned locally with exponential precision. Since neighboring qubits can be coupled and 
decoupled individually, it is sufhcicnt to study and understand the physics of single and 
double quantum dots together with the coupling mechanisms to the environment present 
in particular systems [41]. Undesired interactions between three, four, and more qubits 
should then not pose any great concern. This is in contrast with proposals that make 
use of long-ranged interactions (such as dipolar coupling), where scalability might not 
be easily achieved. 

Figure 2 displays part of a possible implementation of a quantum computer. Dis- 
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played are four qubits represented by the four single electron spins confined vertically 
in the heterostructure quantum well and laterally by voltages applied to the top gates. 
Initialization of the quantum computer could be realized at low temperature T by apply- 
ing an external magnetic field B satisfying |.9/iB-B| ^ kgT, where g is the g-factor, /i^ 
is Bohr's magneton, and ks is the Boltzmann constant. After a sufficiently long time, 
virtually all spins will have equilibrated to their thermodynamic ground state |0) = |t). 
As discussed in the 2nd criterion of the last section, this method might be too slow for ze- 
roing qubits in a running computation. Other proposed techniques include initialization 
through spin-injection from a ferromagnet, as has been performed in bulk semiconduc- 
tors [42, 43], with a spin-polarized current from a spin-filter device [44, 45, 7, 46, 47], or 
by optical pumping [48, 49, 50, 51]. The latter method has allowed the preparation of 
spin states with very high fidelity, in one case as high as 99.8% [52]. 

The proposal of ref. [7] requires single qubit rotations around a fixed axis in order to 
implement the CNOT gate (see below). In the original work [7] this is suggested to be 
accomplished by varying the Zeeman splitting on each dot individually, which was pro- 
posed to be done via a site-selective magnetic field (generated by, e.g., a scanning-probe 
tip) or by controlled hopping of the electron to a nearby auxiliary ferromagnetic dot. 
Local control over the Zeeman energy may also be achieved through g-factor modulation 
[53], the inclusion of magnetic layers [54] (see also fig. 2) or by modification of the local 
Overhauser field due to hyperfine couplings [55] . Arbitrary rotations may be performed 
via ESR induced by an externally applied oscillating magnetic field (sec subscct. 3 '4). 
In this case, however, site-selective tuning of the Zeeman energy is still required in order 
to bring a specific electron in resonance with the external field, while leaving the other 
electrons mitouchcd (see also fig. 2). Alternative all-electrical proposals {i.e., without 
the need for local control over magnetic fields) in the presence of spin-orbit interaction 
or a static magnetic field gradient have been discussed recently. See subsubsect. 3'4.1 
for greater details. 

Two-qubit nearest-neighbor interaction is controlled in the proposal of ref. [7] by 
electrical pulsing of a center gate between the two electrons. If the gate voltage is high, 
the interaction is 'off' since tunneling is suppressed exponentially with the voltage. On 
the other hand, the coupling can be switched 'on' by lowering the central barrier for a 
certain switching time t, . In this configuration, the interaction of the two spins may be 
described in terms of the isotropic Heisenberg Hamiltonian 

(1) Hs{t) = J{t)SL ■ Sr, 

where J{t) oc tQ{t)/U is the time-dependent exchange coupling that is produced by 
turning on and off the tunneling matrix element to{t) via the center gate voltage. U 
denotes the charging energy of a single dot, and and are the spin-i operators for 
the left and right dot, respectively. Equation (1) is a good description of the double-dot 
system if the following criteria are satisfied: (i) AE ksT, where T is the temperature 
and AE the level spacing. This means that the temperature cannot provide sufficient 
energy for transitions to higher-lying orbital states, which can therefore be ignored, (ii) 
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Ts ^ /S.E/h, requiring the switching time Tg to be such that the action of the Hamiltonian 
is 'adiabatic enough' to prevent transitions to higher orbital levels, (iii) U > to{t) for all 
t in order for the Heisenberg approximation to be accurate, (iv) F^^ 3> t^, where 
is the decoherence time. This is basically a restatement of the 3rd DiVincenzo criterion. 
For recent experimental results on the decoherence times in lateral GaAs quantum dots, 
see subsubsect. 3 '3. 3. 

The pulsed Hamiltonian eq. (1) applies a unitary time evolution Us{t) to the state of 
the double dot given by Us{t) = Texp[-i H,{t')dt' /h] = exp[-{i/h) J{t')dt'SL-Sii]. 
If the constant interaction J{t) = Jo is switched on for a time Ts such that J^" J{t)dt/h = 
JoTs/h = TT mod 27r, then Us{ts) exchanges the states of the qubits: Us(Ts)\n,n') = 
|n',n). Here, n and n' denote real unit vectors and |n, n') is a simultaneous cigenstate 
of the two operators Sl • n and S_r • n'. This gate is called SWAP. If the interaction is 
switched on for the shorter time Ts/2, then Us{ts/2) = Us{tsY^'^ performs the so-called 
'square-root of swap' denoted by VswAP. This gate together with single-qubit rotations 
about a fixed (say, the z-) axis can be used to synthesize the CNOT operation [7] 

(2) C/coT = e^(-/2)^^e-'(-/2)^«C/,(T,)V2e-^£c/,(r«)V2, 
or, alternatively, as 

(3) Uc.o. = e-^^C/«(T,)-V2e-'(-/2)5ic/,(r,)e'(-/2)^£c/,(r,)i/^ 

The latter representation has the potential advantage that single qubit rotations involve 
only one spin, in this case the one in the left dot. Writing the CNOT gate as above, it is 
seen that arbitrary single qubit rotations together with the v^SWAP gate arc sufficient for 
universal quantum computing. See subsect. 3'4 for a recent experimental implementation 
of the \/swAP operation. Errors during the execution of a \/swap gate due to non- 
adiabatic transitions to higher orbital states [56, 57], spin-orbit interaction [58, 59, 60], 
and hyperfine coupling to surrounding nuclear spins [61, 62, 63, 64] have been studied. 
Furthermore, realistic systems will include some anisotropic spin terms in the exchange 
interaction which may cause additional errors. Conversely, this fact might be used to 
perform universal quantum computing with two-spin encoded qubits, in the absence of 
single-spin rotations [58, 65, 66, 67]. 

2'6. Alternative approaches to quantum, com,puting . Although the remainder of the 
review will mostly be concerned with the realization of the spin-qubit proposal of ref. [7] 
(and related decoherence effects), we would nevertheless like to discuss some of the al- 
ternative proposals for quantum computing that have emerged in recent years. Note 
that by this we are not referring to the many alternative physical implementations of 
qubits which are also studied extensively in present-day research (see, e.g., ref. [68] 
for a review focussing mainly on solid state qubits). Rather, we would like to review 
proposals for quantum computers which fundamentally differ from the standard circuit 
model. The schemes we will turn our attention to are measurement-based and adiabatic 
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quantum computing. We will not discuss topological quantum computing in greater 
detail, which performs computation by braiding non-Abelian anyons. These are par- 
ticular quasi-particle excitations predicted to exist in certain two-dimensional strongly 
correlated many-body systems such as a two-dimensional electron gas in the fractional 
quantum Hall regime. Topological quantum computing is supposed to be much less sus- 
ceptible to gate errors since small deformations of braids do not change their topology. 
The interested reader is referred to the recent reviews refs. [69, 70]. Sec ref. [71] for a 
measurement based implementation of CNOT on = 5/2 Ising-type anyon qubits. 

2'6.1. Measurement-based quantum computing. Implementing quantum gates, partic- 
ularly two-qubit gates, with a precision as required by fault-tolerant error correction is 
difficult. Instead of performing gate operations on qubits, there are proposals that allow 
for miiversal quantum computing by replacing part or all of these gates by measurement. 
We will mainly focus on a measurement-based implementation of CNOT for qubits repre- 
sented by single or multiple electron spins. Afterwards, we will briefly outline the ideas 
behind the so-called 'one-way quantum computer'. 

Measurement-based implementation of CNOT. When using the polarization state of a pho- 
ton as a qubit, it is known that universal quantum computing can be achieved using only 
linear optics and single photon measurements [72]. This holds similarly for all bosons. 
For electrons (and, similarly, for all fcrmions), there exists a strong no-go theorem [73, 74] 
stating that quantum computing with single-electron Hamiltonians and single-spin mea- 
surements can efficiently be simulated on a classical computer, thus not exhibiting the 
observed exponential speed-up of some algorithms over their classical analogs. How- 
ever, the no-go theorem can be circumvented by exploiting the electron's charge degree 
of freedom: It has been shown recently how to build a CNOT gate for single- [75] and 
multi-electron [71] qubits by the ability to perform, apart from the availability of single- 
electron operations and single-spin readouts, charge measurements. Universal quantum 
computing is thereby restored. Note that the qubits are still encoded in spin states of 
electrons. Since spin and charge are commuting observables, charge measurements do 
not alter the information represented in the spins. 

The main idea is to provide parity measurement of two electron spins via charge 
detection. The cnot gate is then constructed from these parity gates. Imagine we had 
such a device at hand, i.e., for a state in a space either spanned by {|TT))lii)} (even 
parity), or by {jti), jit)} (odd parity), it could determine nondcstructivcly which space 
the state belongs to by detecting the presence or absence of charge upon a measurement 
thereof. A rather abstract notion of a parity gate was described in ref. [75]. Further below, 
we will describe a much more concrete theoretical proposal in the reach of presc;nt-day 
experiments. In the following, however, we will first review the realization of a CNOT gate 
for single-electron qubits as presented in ref. [75]. A generalization of this proposal to 
multi-electron qubits can be found in ref. [71] . Particularly, a detailed construction of the 
CNOT gate for two-electron qubits encoded in the singlet-triplet basis (see subsect. 3'2) 
is discussed there. 
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control in control out |LL) - states |LR) - states |RR) - states 




target in target out jjjj 

ITT) 



Fig. 3. - a) Construction of a deterministic CNOT gate from two parity gates Pi and P2. Each 
parity gate has two input and two output arms and can discriminate whether two incident 
electrons are in a parallel or antiparallel spin configuration via charge sensing. The gates H 
symbolize Hadamard transformations. The CNOT operation is completed by applying Pauli 
rotations to the control and target electron depending on the outcomes of the two parity and 
the ancilla measurements. See main text for a detailed discussion. (Figure adapted from ref. 
[75].) b) Energy diagram for the double-dot parity gate proposed in ref. [76]. The charge states 
]LL), ]LR), and |RR) are indicated by the electron configuration in the double dot. The energy 
levels of all possible spin states are shown for each charge configuration. Gate voltages and a 
magnetic field are applied such that the ]LR) states are lowest in energy, and the antiparallel 
]LL) and |RR) states are on the same energy. This allows for resonant tunneling between the 
latter two states (indicated by the dotted arrow), while tunneling between parallel states is 
suppressed due to the Zeeman mismatch A§ > A^. (Figure adapted from ref. [76].) 



Let a parity gate work as follows. Two electrons can enter the gate simultaneously, 
and after the parity was measured via charge detection, the electrons leave the gate with 
unmodified spin state if the latter was in one of the even or odd parity spaces described 
above. Furthermore, let the gate record 'no charge' {p = 0) for two antiparallel spins, 
and 'charge' {p — 1) for parallel incident spins. Figure 3 a displays the construction of 
the deterministic CNOT gate using two connected parity gates Pi and P2- Before the 
input and after the output arms of P2, a Hadamard transformation H is applied to each 
spin, defined as |T) ^ (|T> + \l))/V2, and |i) (|t) - \l))/V2. The control qubit enters 
the first gate Pi. Its state decides whether the target qubit, entering the second gate P2, 
is to be flipped according to the definition of the CNOT operation. Pi is also provided 
with an ancilla qubit prepared in the state (|T) + li))/\/2 which is then fed back into 
P2- Upon leaving the second gate, the ancilla is measured. Conditioned on this result 
and the outcomes pi and p2 of the two parity measurements in Pi and P2 , respectively, 
a Pauli matrix has to be applied to control and target qubit in order to complete the 
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CNOT operation (see below). We will now see (following the supplementary appendix of 
ref. [75]) that this setup indeed implements cnot. 

Wc depart from our usual notation for the spin basis and identify |0) = |t) and 
|1) = II). Furthermore, all variables represent a number in {0,1} and addition is per- 
formed modulo 2. We first consider the action of the second gate P^- After applying 
the Hadamard gates on the input arms of P2, but before the parity measurement, an 
input state \a)\y) has been transformed to (|0) + (— 1)°|1))(|0) + (— 1)^|1)) (normalization 
constants will be neglected for the rest of this section). Here, the first (second) state 
represents the qubit entering the upper (lower) arm of the parity gate. After the parity 
measurement the state has become 



(4) \a)\y) 



f|0)|0) + (-l)(''+«)|l)|l) if P2 = l 
\(-l)^|0)|l) + (-l)«|l)|0) if P2 = 0. 



In thc> end, the Hadamard gates on the output arms are performed and the state of the 
ancilla is measured, i.e., 



(5) \a)\y) 



|0)|a + 2/) + |l)|a + y + l) if P2 = I 

_(-l)'^|0)|a + 2/)-(-l)''|l)|a + 2/+l) if P2 = 
(6) ^ (-l)(f=+^)('*+^)|o + y + ^), 



where z is the outcome of the ancilla measurement. The action of the first parity gate 
Pi on the control and the ancilla qubit is given by |a;)(|0) + |1)) \x)\x +pi + 1). The 
second state is transmitted to the upper arm of P2, yielding the total action of the setup 
on a control-target pair Is)]?/): 

(7) \x)\y) ^ (-l)(f^+i)(^+^+fi+i)|x)|a; + y + z+pi + l). 

Post-correction depending on pi, p2 and z has to be performed now in order to ob- 
tain the correct CNOT operation defined by \x)\y) \x)\x + y). The phase factor 
(_l^(j'2+i)(z-i-j)i+i) jg irrelevant since it does not depend on x and y. If p2 = 0, a az 
gate has to be applied to the control qubit. This eliminates the remaining phase (since 
az\x) = (— l)^|a;)). In order to obtain the correct target \x + y), a cfx gate needs to be 
applied li z+pi =0 (since (Jx\y) = |y + !))• This completes the description of the CNOT 
gate in terms of parity measurements. 

We now qualitatively describe a concrete proposal due to ref. [76] of a parity gate 
exploiting charge measurement. The device consists of two coupled quantum dots con- 
taining the two electrons whose parity is to be determined. The dots are assumed to 
have different Zeeman splittings and A^. This could be realized, e.g., by locally 
different magnetic fields or with an inhomogeneous g- factor. By applying suitable gate 
voltages and a perpendicular magnetic field, one can achieve the energy configuration 
shown in fig. 3b. The gate voltages are set such that all states |LL) (two electrons in 
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the left dot) and |RR) (two electrons in the right dot) are higher in energy than |LR) 
(one electron in each dot), independent of the spin configuration. The strength of the 
external magnetic field is chosen such that the zero-field singlet-triplet splitting in each 
dot is removed (see subsubsect. 3'2.1). This leads to the degeneracy of all spin states 
in the odd-parity space {|Ti))|iT)} with charge configuration |RR) or |LL). The gate 
voltages can further be properly tuned to align these degenerate levels of the two dots. 
However, due to the different Zceman splittings in the left and the right dot, parallel 
spin configurations with charge state |RR) are detuned by e = ±2(A§ — A^) from the 
corresponding states in |LL). The energy spectrum hereby achieved allows for elastic 
tunneling between the states |LL) and |RR) (through the intermediate state |LR)) with 
antiparallel spins, whereas resonant tunneling for parallel spins is suppressed due to the 
Zeeman mismatch. The transition to the ground state occurs only inelastically. A quan- 
tum point contact near the neighboring dots is used as an electrometer to detect the 
presence (indicating antiparallel spins, as opposed to the abstract gate described above) 
or absence (indicating parallel spins) of tunneling events (see also subsect. 3 '5). A mi- 
croscopic model of the double dot system is further studied in ref . [76] , where it is shown 
that elastic tunneling (if present) strongly dominates over inelastic tunneling, and that 
the device still works with high fidelity even if the measurement parameters cannot be 
controlled perfectly. 

The one-way quantum computer. A proposal for quantum computing that requires noth- 
ing but single qubit measurements during computation is briefly outlined in the following. 
The one-way quantum computer [77] requires a so-called cluster state [78] to begin with. 
This is a certain highly entangled state that can be realized on a two- or three-dimensional 
array of qubits interacting through externally controllable nearest-neighbor Ising- [78] or 
Heisenberg-type [79] interactions. After this state is initialized, a network, i.e., an entan- 
gled state among qubits forming a grid-like structure, is realized by discarding undesired 
qubits through measurements of az. The computation is then performed by measure- 
ments in the x-y plane of qubits in the network. The choice of future measurement bases 
may depend on past measurement outcomes. It can be shown that universal quantum 
computing is possible, whereby the computation proceeds spatially from left to right 
with quantum information flowing on horizontal branches on the network and two-bit 
interactions implemented by measurements on vertical branches. The interested reader 
is referred to refs. [77, 78, 80] for greater detail. 

Small-scale instances of Grover's algorithm [81, 82] and the Deutsch-Jozsa algorithm 
[83] have recently been demonstrated within quantum optics. 

2"6.2. Adiabatic quantum computing. Adiabatic quantum computing appeared flrst 
in the context of a novel approach to solve classical optimization problems [84]. It 
has then evolved into a general approach to quantum computation now known to be 
polynomially equivalent to the standard model [85] , implying that standard and adiabatic 

quantum computers have the same computational power. While it was found rather 
quickly that a standard quantum computer can efficiently simulate arbitrary adiabatic 
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computations [86, 87], hence proving one direction of the equivalence, it took several 
years to show the opposite, i.e., that an adiabatic quantum computer can simulate any 
standard computation with only a polynomial overhead [85]. We will review the basic 
ideas behind adiabatic quantum computing and the original approach towards it [85]. 
Afterwards, we will point out some quite recent developments in the field. 

The basis of adiabatic quantum computing is the adiabatic theorem [88, 89] : Given a 
system initially prepared in an energy eigenstate and undergoing an externally induced 
time evolution, the theorem says that the system's state will remain arbitrarily close to 
the corresponding instantaneous eigenstate, if there is a nonzero energy gap all along the 
evolution and if the latter is carried out 'slow enough'. Hereby, the time scale for 'slow 
enough' depends on the desired closeness accuracy of the system's state to the respective 
instantaneous eigenstate and on the size of the minimal gap along the evolution. The 
smaller the minimal gap, the slower the process has to be performed in order to suppress 
transitions to states higher in energy. In adiabatic quantum computing, one starts with 
a system prepared in the ground state of an initial Hamiltonian -ffinit, with that ground 
state being unique and having a simple form such as jOOO . . .). The system is then being 
evolved adiabatically according to H{t) = (1 - t/T)Hinit + {t/T)Hfined, t e [0, T]{^), into 
a setup described by the final Hamiltonian -ff final, whose ground .state encodes the result 
of the desired computation. This has the potential advantage of not requiring fast gate 
and measurement operations. Additionally, adiabatic quantum computing is intrinsically 
robust against decoherence due to environmental noise (see ref. [91] for a review). One 
can prove a rigorous lower bound on the value of T required to obtain a final state that is 
£-close in i2-norm to the ground state of iffinai- T depends inversely both on e and on the 
minimal gap between the ground and first excited state of H{t) (called the spectral gap). 
The running time r is defined to be Tmaxt ]]ff(t)]], where the second factor makes r 
invariant to rescaling ofH{t). Further, ffinit and ifsnai are restricted to be, in the original 
terminology, local, meaning that they may only allow interactions between a constant 
number of particles in order to be physically realistic. This constraint also assures that 
the Hamiltonians have efficient classical descriptions [85]. 

The main difficulty to adiabatic quantum computing is the fact that finding an iJfinai 
encoding the result of a computation in its ground state is impossible, since that result is 
intrinsically unknown (otherwise there would be no need for its computation). While in 
the standard model an algorithm is executed by a discrete unitary time evolution, iJfi„ai 
is subjected to simultaneous local constraints. This problem is overcome in ref. [85] by 
loosening the requirement for //final to have a ground state exactly equal to the outcome 
of a standard quantum computation. It is shown that it is sufficient to obtain a ground 
state having nonzero overlap with the desired state. This overlap can then be enhanced 
arbitrarily with polynomial overhead. 

Given a quantum circuit with L gates and denoting by \a(l)), I = 0,1,..., L, the 
state in the circuit after the Zth gate Ui has been applied, it is aimed at constructing a 



(^) General non-linear paths have been studied as well, see, e.g., ref. [90]. 
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Hamiltonian ifgnai whose ground state is the so-cahed history state 



1 



L 



(8) 



The right L qubits are referred to as clock qubits (superscript 'c') whose representation 
"enables a local verification of correct propagation of the computation from one step 
to the next, which cannot be done without the intermediate computational steps" [85]. 
After showing that there is an iJinit with non-degenerate ground state |00 . . . 0) ® 
-fffinai is defined as 



where -ffjnput and ffciock ensure that undesired input states and illegal clock states receive 
an energy penalty. The Hi, 1 < I < L are given by 



and similar for / = 1 and I = L. The Hi make sure that the unitary action of each 
Ui comes along with the correct update of the clock register. The final Hamiltonian 
constructed in this way has indeed jry) as its ground state. In ref. [85] it is then further 
shown that the spectral gap of H{t) is lower bounded by an inverse polynomial in L 
for all t. This, together with adding identity gates (in the circuit picture) at the end of 
the computation in order to arbitrarily increase the weight of \oi{L)) in \r]), results in a 
running time that scales with L^. 

Note, however, that the Hamiltonian obtained in this way is 5-local, i.e., interactions 
between 5 arbitrarily distant particles have to be realized. It was already shown in 
ref. [85] that -fffinai can be made 3-local, but the running time then increases and roughly 
scales as L^^. It was even demonstrated how to make an iJgnai involving only two- 
body nearest-neighbor interactions, although this required the usage of particles having 
a six-dimensional state space. These results have recently been extended to qubits with 
2-local interactions [36, 92] and to qubits on a 2D lattice with nearest-neighbor two-body 
interactions [93]. All these constructions start from the 5-local Hamiltonian described 
earlier. Very recently, a rather different approach to adiabatic quantum computing has 
been taken using the concept of 'ground state quantum computation' [94]. There, the 
entire temporal trajectory of an algorithm is encoded spatially in the ground state of a 
suitable system. The proposal uses qubits and requires only two-body nearest-neighbor 
interactions. It has been shown that the scaling of the running time with N and L, where 
A'^ is the number of qubits, is of order (NL)'^. The interested reader is referred to the 
original literature [94, 95]. 



(9) 




1=1 



(10) 



Hi=I^ |100)(100|f_i,,,;+i -Ui® |110)(100|F_i,i,,+i 
- U} ® |100)(110|l^_i,;,;+i +7® |110)(110|ti,;,;+i, 
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2"7. Entanglement measures. - Quantum correlations are heavily exploited in every 
quantum algorithm and form the key ingredient to the reason why quantum computing 
differs from classical computation. In the circuit model discussed in subsect. 2"2, en- 
tanglement is generated by the CNOT gate. Typical simple examples of entangled states 
are Bell states such as |$"'") = (|t)A|T)B + |i)A|i)B)/'\/2 [96]. Imagining that two par- 
ticles in such a state travel to spatially arbitrarily separated observers A and B, both 
observers will obtain random, but perfectly correlated measurement outcomes (assum- 
ing ideal measurement conditions). This is inconsistent with any classical {i.e., local) 
description of the state [96, 97] . 

Entanglement hence manifests itself in the form of inter-partite correlations in a quan- 
tum state which arc not explainable by classical means. It this context, one usually 
introduces the notion of LOCC-operations [98]: If correlations observed in a quantum 
state cannot be reproduced (or simulated), starting from initially unrelated quantum 
subsystems, by local quantum operations ('LO') coordinated by and influencing each 
other via classical communication ('CC'), then these correlations are identified with the 
presence of entanglement in that state. On the other hand, if a state can be created by 
LOCC-operations alone, it is denoted separable, i.e., unentangled. The fact that LOCC- 
operations can neither create entanglement in a separable state, nor enhance already 
present entanglement (on average(^)), makes entanglement a resource which is sought to 
be quantified. 

Formally, a pure state in an n-partite Hilbert space H = Hi is called entangled 

if it cannot be written as a product state \tp) = ^^^i {fpi). For example, the state = 
(ITT) - ITi) + ITT) - Ui))/2 can be written in the form |0) - (|t) + |i))/V2® (|t) - |i))/V2 
and is thus separable, i.e., not entangled, whereas such a decomposition is not possible 
for the Bell state |$+) from above. Analogously, a mixed state p acting on H is separable 
if it can be written in the form p = "^^Pi (S'fcLi where the p\. act on Tik for all k. 
However, as already indicated earlier, the story is not over after categorizing states into 
a 'black and white' scheme by determining whether a particular state is separable or 
entangled. A simple example demonstrating that some states can be 'more entangled' 
than others is the fact that some states violate Bcll-type inequalities stronger than others, 
implying the presence of more quantum correlations. 

An entanglement measure is a function from the space of density matrices to a closed 
interval in the real non-negative numbers, the lower bound usually being 0, and should 
reflect the physical properties of entanglement. Most importantly, it should be non- 
increasing under LOCC-operations on average, which is meaningful, taking the previous 
discussion into account. In particular, it should be invariant under local unitary trans- 
formations which merely correspond to local changes of basis. An entanglement measure 



(®) There is a protocol known as the 'Procrustean method' or 'entanglement gambling' [99, 100]: 
One can turn any multipartite entangled state into a Bell state shared by some pair of parties 
using LOCC alone. This works however only with probability smaller than 1. On average, 
LOCC cannot increase entanglement. 
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satisfying the previous conditions is called an entanglement monotone [101, 102]. Addi- 
tionally, entanglement measures are often demanded to be able to uniquely distinguish 
between separable and entangled states, usually incorporated by constructing the mea- 
sures such that they are if and only if the state examined is separable. 

Entanglement in bipartite systems is the case understood by far the most until now, in 
contrast to multipartite entanglement (see, e.g., ref. [98]). This is also due to the existence 
of a meaningful entanglement measure for bipartite states, namely, the entanglement of 
formation Ep [102], defined as 

(11) Ef{p)= min 

{Pi,l^i>}ez)(p)"^ 



where 



(12) 



S)(p) ={{pi, \K > rankp, ft > 0, 

K K 

'Y^Pi = 1, IV'i) e n, /9 = ^Pi|V'i)(V'i|} 



i=l 

is the set of all so-called pure-state decompositions of p, and 

(13) E{m = - tr [(tn |^)(V'|) log2(tri |^)(V'|)] 

is the entropy of entanglement (tri denotes the partial trace over the first subsystem). 
The latter is an entanglement monotone for bipartite pure states and is closely related to 
the von Neumann entropy. The numerical value of Ep{p) is meaningful in the following 
sense: It has been shown [102] that, given a number of N identical states p, one can 
(asymptotically) 'distill' N ■ Ep{p) maximally entangled Bell states (such as \^~^)) out 
of them. The entanglement of formation thus measures quantum correlations in units of 
the entanglement contained in a Bell state(^). 

There exists a vast amount of proposed entanglement measures for multipartite pure 
states. The study of mixed-state entanglement is, however, important as well, since any 
realistic quantum system will eventually couple to the environment and thus decohere. 
The so-called convex roof construction [103] (the entanglement of formation being an 
early example thereof) gives a general recipe how to extend a pure-state entanglement 
monotone to mixed states: Given an arbitrary pure-state multipartite entanglement 
monotone m, the convex roof of m is given by 

(14) M{p)= , inf VftmdV.)), 



(^) Note that for states in highcr-dimcnsional systems, the entanglement of formation can be 
larger than 1, implying that more than one Bell state is required to create such a state. 
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where D{p) is defined as in eq. (12). M has the desirable feature that it is an entan- 
glement monotone itself, and that it properly reduces to m if p describes a pure state 
[104]. The optimization problem coming along with eq. (14) is, however, rather involved 
and seems impossible to bo analytically solvable in general. Remarkably, there is one 
major exception to this statement: There exists a general analytical expression for the 
entanglement of formation of two qubits [105] . 

Nevertheless, the optimization problem in cq. (14) can be tackled numerically to some 
extent [106, 107, 108] by first parameterizing the set of all pure-state decompositions 
D{p). Let St{k,r) denote the set of all fc x r matrices U G c*^^'' with the property 
WU = Irxr- The required parametrization of T){p) is due to the Schrodinger-HJW 
theorem [109, 110], stating that every decomposition of p into k states is related to a 
matrix U € St{k, r), where r = rankp, and vice versa. A search over St{k, r) for all /c > r 
is thus equivalent to searching over S)(p). Explicitly, given a matrix U G St{k, r) and the 
eigendecomposition p = J2^=i ^i\Xi){Xi \ of P) the pure-state decomposition corresponding 
to U is given by 

(15) Pi = mi>i), m = {i/vm)m, 



where 

r 

(16) m = J2^ijV>^\xj), i = i,...,k. 

One is hence confronted with the new optimization problem 

(17) M{p)=min mf hiU,p), 

k>r UeSt{k,r) 

where h{U,p) is the convex sum on the right-hand side of eq. (14). In practice, one 
can of course only investigate this problem for a few values of k. However, numerical 
studies show that a k not much larger than r is sufficient for obtaining accurate results 
[106, 107, 108]. General-purpose numerical algorithms tackling eq. (17) for arbitrary 
pure-state entanglement monotones have been presented very recently [107, 108], together 
with studies of entanglement in states emerging from physical Hamiltonians. 



3. — Spin manipulation in GaAs quantum dots 

As announced, we will restrict ourselves in the rest of the review mostly to lateral 
GaAs quantum dots. This is motivated by the remarkable latest achievements we have 
witnessed in the field. The main features of the spin-qubit proposal of ref. [7], in particular 
the single and two-qubits gates, have by now been realized in single and double GaAs 
quantum dots with various degrees of accuracy. It is therefore meaningful to review the 
field in the light of the five DiVincenzo criteria for scalable quantum computing presented 
in subsect. 2 '4. 
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3'1. Realization of well-defined spin qubits. - The qubits we are considering are ob- 
tained in a standard two-dimensional electron gas (2DEG) formed at the interface of a 
GaAs/AlGaAs heterostructure, as illustrated in fig. 2. The electron gas is then depleted 
by means of metallic top gates in order to define the confinement region of the quantum 
dots. We refer to refs. [40] and [41] for a detailed discussion of the stability diagram 
of single and double quantum dots. The main features thereof can be recovered from a 
simple charging Hamiltonian [41] 

(18) Ho = ^Yl - 1) + ^12^1^2 -eJ2^i^i + Yl ^rnUi^m, 

i i i,m 

where we assumed the two dots (i = 1, 2) to be identical for simplicity. Here, U and U12 
are the on-site and nearest-neighbor Coulomb repulsions, respectively, Vi are the local 
potentials at each dot, and the denote single-particle orbital energies with occupation 
ni^m = ^^j.m.T + ^i,m.j.- Furthermore, Ni = X^„j JT-i.m is the number of electrons in dot 
i. Smaller but still relevant corrections to Hq are determined by external fields {e.g., 
magnetic fields), spin-couplings to the environment {e.g., via spin-orbit or hyperfine 
interaction, sec subscct. 3'3 and sect. 4) and tunneling. The simplest example of such 
an additional term is the Zeeman coupling 

(19) Hz = g,iBB-J2^j, 

3 

where B is the externally applied magnetic field and Sj = (jjjl is the spin-^ operator of 

the j-th electron in the double dot, j = 1, 2, . . . , TVi + 7V2- In the following, we generally 
define the quantization direction to be along B, with the |t) orientation having lower 
Zeeman energy since 5 < in GaAs. 

The metallic gates allow one to control the potentials Vi, which determine the ground- 
state occupation denoted by (A^i, N^)- The spin qubits arc realized by occupying each dot 
with exactly one electron. While control of the electron number down to single occupancy 
was achieved early on for other types of dots {e.g., vertical quantum dots [111]), the lateral 
confinement tends to suppress the tunneling rates with the reservoirs. This problem leads 
to difficulties in observing the few-electron regime but can be overcome by designing 
proper gating structures (figure 4 shows two examples of actual samples). For this reason, 
the first demonstrations of few-electron single [112] and double dots [113, 114, 115] with 
lateral gating are much more recent than for vertical dots. 

3'2. Initialization of the spin state. - A straightforward procedure to initialize a sys- 
tems of qubits is to apply a sufficiently large magnetic field and wait for relaxation to 
the ground state jtTT • • •) to occur. Experiments are usually performed in dilution re- 
frigerators with base temperature around 20 mK, which is smaller than typical Zeeman 

splittings (^ 300 mK at i? = 1 T and using the bulk value g = —0.44). The initialization 
time is of the order of a few relaxation times, which in GaAs dots have been reported 
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200 nm 200 nm 

Fig. 4. - Scanning electron micrographs of a single (left) and double (right) quantum dot. The 
electrons are confined to the regions indicated by the white circles, surrounded by metallic top 
gates (light gray areas), and the square boxes indicate the Ohmic contacts. In the double dot 
setup, two quantum point contacts used for charge sensing can be seen. (Reprinted figures with 
permission from ref. [40]. Copyright (2007) by the American Physical Society.) 



to be as high as ^ 1 s (see subsubsect. 3'3.3 for a more complete discussion). In the 
following, we discuss several other techniques used in practice to initialize single and 

double GaAs quantum dots to configurations other than This allows, e.g., for more 
flexibility in the subsequent manipulation of the double dot spin state. 

3'2.1. Singlet-triplet transition in single dots. We consider here an isolated dot (more 
precisely, the i = 2 dot in eq. (18)) and show that, if two electrons are present, the 
ground state can be chosen to be either a singlet or a triplet, depending on the value of 
the external magnetic field. Initialization in the desired spin state can, in principle, be 
accomplished easily by energy relaxation. We start from the lowest energy single electron 
states with charge configuration (0,1). Clearly, they only differ due to the spin and have 
energies £'±(0, 1) = eo - eV-i T (see cqs. (18) and (19)), where l^^Ez = \gpiBB\ is 

the Zeeman splitting. If now one more electron is added, singlet and triplet states can be 
formed. The lowest lying states are denoted by 5(0, 2) and Ti„(0, 2), where the subscript 
m = 0, +, — refers to the component of the total spin parallel to B. The energies are 
given by £<?(0, 2) = 2(eo - eV2) + U and £'t,„ (0, 2) = Es{Q, 2) + I^Est - mAEz, where 
the single-triplet splitting at zero magnetic field is given by AEst = Acorb — . Here, 
^Eorb = ei — Co is the difference in orbital energies, which would be the only contribution 
to AEst according to the simple charging Hamiltonian cq. (18). However, the splitting 
is experimentally found to be smaller than Atorb due to a change AU in the charging 
energy [40]. The splitting AEst is generally still positive, resulting in a singlet ground 
state at zero magnetic field. 

Interestingly, the Zeeman energy is smaller than AEst for typical values of the mag- 
netic field. Therefore, an in-plane magnetic field cannot induce a singlet-triplet transition 
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Fig. 5. - In the left panel, the Pauli spin blockade is illustrated. At positive bias, an electron 
from the left reservoir can tunnel to the right contact only if it enters the (0,1) double dot in a 
singlet configuration 5(1, 1), while the triplet ^(l, 1) is blocked due to spin conservation (dashed 
transition). At negative bieis the electron can only tunnel in the S{0, 2) singlet and transport is 
always possible. The position of the chemical potentials /ia(l, 1) (a ~ S,Tm) is indicated in the 
first dot, while /xs(0, 2) and /xt„(0, 2) are shown in the second one. The detailed structure of 
the (1,1) and T(0, 2) levels is not specified, since it is of an energy scale much smaller than the 
5(0, 2)-r(0, 2) singlet-triplet splitting (see subsubsect. 3'2.3 for a more detailed discussion). 



since, as a good approximation, it does not affect the orbital states. Instead, a significant 
decrease of Acorb and increase of MJ is produced by a magnetic field B±_ perpendicular 
to the 2DEG due to orbital effects. The energy crossing of singlet and triplet is typically 
realized around ~ 1 T [116, 117]. This condition is required for the realization of 
the parity gate discussed in subsubsect. 2'6.1. Furthermore, 5° can be tuned via elec- 
tric gates [117]. In the following, however, we will usually neglect orbital effects of the 
magnetic field, supposing B to be either in-plane or sufficiently small. 

3'2.2. Pauli spin blockade in double dots. The lowest-lying (1,1) spin states are the 

singlet S'(l, 1) and the triplets Tm(l, 1). The energies i?Q(l, 1) (where a = Tm) are 
degenerate in first approximation and eq. (18) gives i?a(l) 1) — E{1, 1) = 2eo— e(Vi-|-V2)-|- 
U\2- It is however possible to selectively prepare the system in a triplet state via Pauli spin 
blockade [118, 119]. This is realized at positive bias if the chemical potentials /ia(l) 1) and 
/Xa(0, 2) of the double dot are adjusted as shown in the left panel of fig. 5. The chemical 
potentials with respect to the (0, 1) occupation are defined as ^a(0, 2) = Ea{Q, 2)~-E{Q, 1) 
and /ic«(lil) = -E'q(1,1) — £^(0,1), and are shown in fig. 5 (for simphcity, we neglect 
the presence of small Zeeman splittings and assume £^+(0,1) ~ i?_(0, 1) = E{{),1)). 
Tunneling occurs from the left reservoir which is connected to the first dot, to the right 
reservoir connected to the second dot. The sequence (0, 1) — > (1, 1) — > 5(0, 2) — > (0, 1) 
would be energetically allowed both through ^(l,!) and T„j(l,l), but the transition 
Tto(1,1) — > 5(0,2) is forbidden due to spin conservation. Therefore, as soon as an 
electron tunnels from the left reservoir to ^^(l, 1), the double dot is blocked in the 
triplet state. Transport is only possible after relaxation into 5(1, 1), which can occur on 
a millisecond time scale. 

Note that at negative bias (cf. the right panel of fig. 5) a finite current can flow 
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Fig. 6. - Double-dot energy diagram as a function of the detuning e, close to e = 0. Without 
tunneling (left panel) 5(1, 1) and ro(l, 1) are always degenerate, while the T±{1, 1) are split by 
a finite magnetic field. At e = 0, the two singlet states S{1, 1) and S(0, 2) are also degenerate. 
Inter-dot tunneling (right panel) causes mixing of the ^(l, 1) and S{0, 2) states. The effect of a 
small difference A_Bjv in the nuclear fields of the two dots is not shown here, and is only relevant 
when the splitting between S{1, 1) and Tb(l, 1) is small (i.e., in the left panel (independent of 
e) and at large detuning (je| ^ 0) in the second one). Then, the appropriate eigenstates have 
spin Itl) and lit), as discussed in the text. 



through the dots following the sequence (0, 1) — > S{0, 2) — > S{1, 1) — > (0, 1). In this case 
electrons can only tunnel from the right reservoir to a singlet state, since the triplet is 
too high in energy, and the Tm(l, 1) states are never involved. The Pauli spin blockade 
effect thus leads to current rectification. 

3 '2. 3. Singlet-triplet and charge states mixing in double dots. In the presence of an 
external magnetic field, the triplet energies are i?T,„ — E{1, 1) — mAEz- Therefore, the 
states To(l, 1) and S{1, 1) remain degenerate. However, the magnetic fields on the two 
dots generally differ slightly due to different nuclear configurations (see subsubsect. 3 '3. 2 
on the hyperfine interaction). We assume that the value of the magnetic field is BzLAB^ 
at the left/right dot, where ABn ^ 2 mT. This causes the relevant eigenstates to 
have spin configuration or jjl) (singlet-triplet mixing), with energies = 
E{1,1) =F gfiBAB]^. A first consequence of this fact is that the Pauli spin blockade 
discussed above occurs only for a mixture of T±(l, 1), but not for To(l, 1). The reason is 
that To{l,l) rotates to S'(l, 1) due to the field inhomogeneity ABn, which removes the 
spin blockade since 5(1, 1) can tunnel to S{0, 2). A second consequence is the possibility 
to initialize the system into the spin configurations |ti) or |J,t)- This requires a more 
sophisticated procedure relying on mixing of charge states, as described in the following. 

For simplicity, we neglect for the moment the small field inhomogeneity due to AB^. 
We also define the detuning e of the local potentials at the two dots as follows: eVi = 
eVi — e/2 and eV2 = eVz + e/2, where the Vi are constant potentials such that S{1, 1) 
and S{0, 2) are degenerate at e = 0. The detuning e changes the energy of the (0,2) 
singlet, Es{0,2) — Es{l, 1) — e, while Es{l, 1) and other (1,1) states are not affected. 
These double dot levels are shown in the left panel of fig. 6, while the triplet states 
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Tm{0, 2) have much higher energy and are thus not depicted. Consider next the effect of 
tunnehng in the vicinity of £ = 0. The tunnehng Hamiltonian Ht has a matrix element 
{S{1,1)\Ht\S{0,2)) = \pltc and causes mixing of the singlets with different charging 
configurations. At e = perfect mixing is realized, with energy splitting 2\/2tc, while at 
large e the unperturbed eigenstates 5'(0, 2) and 5'(1,1) are recovered. At large positive 
detuning, 5(0,2) is lower in energy than the T+(l, 1) state, and initialization can hence 
be performed via energy relaxation. If one then slowly changes e toward negative values, 
the system evolves adiabatically along the lower singlet branch into the <S'(1, 1) state (cf. 
the right panel of fig. 6). The leakage to the 5(0, 2) state is estimated in ref. [41]. Note 
also that ^^(l, 1) is not mixed with the singlet in this simple model(^). 

We now consider the effect of ABjv, which is important whenever the S — Tq splitting 
is small. Since it usually holds that tc ^ l^/zsASivl, this is only the case if the detuning 
becomes large in magnitude. In this limit the splitting goes to zero (see fig. 6) and the 
inhomogcneity becomes the dominant effect. If the detuning e is decreased from large 
positive to large negative values faster than the time scale determined by ABn, the 
system will be initiahzed to <S'(1, 1) and will then begin to oscillate between <S'(1, 1) and 
To(l, 1) with frequency 2giiBABN /h. Instead, by adiabatically reducing the value of e, 
the system can be initialized to the spin configuration with lower energy, W\) or 
depending on the sign of ASjv. 

3'3. Relaxation and decoherence in GaAs dots. - The requirement of sufficiently long 
coherence times is perhaps the most challenging aspect for quantum computing architec- 
tures in the solid state. It requires a detailed understanding of the different mechanisms 
that couple the electron's spin to its environment. We introduce here the main concepts 
relevant for GaAs dots, while a detailed discussion is postponed to sect. 4. 

3'3.1. Spin-orbit coupling. While fluctuations in the electrical environment do not 
directly couple to the electron spin, they become relevant for spin decoherence in the 
presence of spin-orbit interaction. In GaAs 2DEGs two types of spin-orbit coupling 
(Dresselhaus and Rashba) are present. The Dresselhaus spin-orbit coupling originates 
from the bulk properties of GaAs [120]. The zinc-blend crystal structure has no center 
of inversion symmetry and a term of the type Hff cx Px{Py — pD'^x + Pyipl ^ Px)'^y + 
Pz{Px ~ Py)'^z allowed in three dimensions, where p is the momentum operator and 
cr are the Pauli matrices. Due to the confining potential along the ^-direction, we can 
substitute the Pz operators with their expectation values. Using (p^) ^ and (pz) = 0, 
one obtains 

(20) Hd = P{Py(^y - Pxl^x)- 



(*) In reality, small spin perturbations cause anticrossing of the singlet branch with T+(l, 1). 
In experiment, e is swept faster around the S — T+ degeneracy in order to avoid the T+(l, 1) 
state [61]. 
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Smaller terms cubic in p have been neglected, which is justified by the presence of strong 
confinement. 

The Rashba spin-orbit coupling is due to the asymmetry of the confining potential 

[121] and can be written in the suggestive form Hr oc {£ x p) ■ a, where £ = £z is an 
effective electric field along the confining direction: 



The Rashba and Dresselhaus terms produce an internal magnetic field linear in the 
electron momentum defined by B so = —2[{Ppx + apy)ex — {Ppy+apx)ey]/giJ,B- If/3 = 0, 
the magnitude of Bgo is isotropic in p and the direction is always perpendicular to the 
velocity. While moving with momentum p, the spin precesses around B<jo and a full 
rotation is completed over a distance of order Ago = \h/{am*)\ = 1—10 ,um, where 
m* is the effective mass. Generally, Rashba and Dresselhaus spin-orbit coupling coexist, 
their relative strength being determined by the confining potential. This results in the 
anisotropy of the spin-orbit coupling in the 2DEG plane (e.g., of the spin splitting as 
function of p). In this case, two distinct spin-orbit lengths can be introduced 



For GaAs quantum dots, the spin-orbit interaction is usually a small correction that 
can be treated perturbatively since the size of the dot (typically ~ 100 nm) is much 
smaller than the spin-orbit coupling lengths A±. The qualitative effect introduced by 
the spin-orbit coupling is a small mixing of the spin eigenstates. As a consequence, the 
perturbed spin eigenstates can be coupled by purely orbital perturbation even if the 
unperturbed states have orthogonal spin components. Relevant charge fluctuations are 
produced by lattice phonons, surrounding gates, electron-hole pair excitations, etc. with 
the phonon bath playing a particularly important role (see subsect. 4"1). 

3 '3. 2. Hyperfine interaction . The other mechanism for spin relaxation and decoher- 
ence that has proved to be effective in GaAs dots, and ultimately constitutes the most 
serious limitation of such systems, is due to the nuclear spins bath. All three nuclear 
species ®^Ga, ^^Ga, and ^'"'As of the host material have spin 3/2 and interact with the 
electron spin via the Fermi contact hyperfine interaction 



(21) 



Hr = a{Px<Jy -PyCTx). 



(22) 



n 



m*{(3 ± a) ' 



(23) 




where Ai and I, are the coupling strengths and the nuclear spin operator at site i, 

respectively. The density of nuclei is no = 45.6 nm~^ and there are typically N ^ 10® 
nuclei in a dot. The strength of the coupling is proportional to the electron density at 
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site i, and one has Ai = A\^p{ri)\'^ /no, where tp{r) is the orbital envelope wave function 
of the electron and ^ w 90 /xeV(^). 

The study of the hyperfine interaction (23) represents an intricate problem involving 
subtle quantum many-body correlations in the nuclear bath and entangled dynamical 
evolution of the electron's spin and nuclear degrees of freedom. While these topics 
will be discussed much more deeply in subsect. 4'2, it is nevertheless useful to present 
here a qualitative picture based on the expectation value of the Ovcrhauscr field Bjv = 
Y^^Aili/gjiB- This field represents a source of uncertainty for the electron dynamics, 
since the precise value of Bjv is not known. Due to the fact that the nuclear spin bath 
is in general a complicated mixture of different nuclear states (see subsubsect. 4"2.1 
for a more detailed discussion of the nuclear density matrix), the operator B^r in the 
direction of the external field B does not correspond to a well-defined eigenstate, but 
results in a statistical ensemble of values. These fluctuations have an amplitude of order 
BN,max/VN ~ 5 mT since the maximum value of Bjv (with fully polarized nuclear bath) 
is about 5 T. 

Finally, even if it were possible to prepare the nuclei in a specific configuration {e.g., 
ITTit ••■))> the nuclear state would still evolve in time to a statistical ensemble on a 
time scale t„tic- Although direct internuclear interactions are present (e.g., magnetic 
dipole-dipole interactions between nuclei), the most important contribution to the bath's 
time evolution is in fact due to the hyperfine coupling itself, causing the back action of 
the electron spin on the nuclear bath. Estimates of the nuclear bath timcscale lead to 
inuc = 10 — 100 jis or longer at higher values of the external magnetic field B [40]. 

3'3.3. Relevant time scales. We provide here a summary of the relevant time scales for 
spin decoherence in GaAs dots. In the Bloch phenomenological description of the time 
evolution, the spin density matrix p = (1 + P • cr)/2 (where P is the spin polarization) 
satisfies 

(24) P = ^B X P-r(P-Po), 

where the tensor Tij is diagonal in a reference frame with the z-axis along B. With this 
choice, the equilibrium polarization is Pq = Pq^z- The time Ti = is the longitudinal 
spin decay time, or spin-flip time, and describes the energy relaxation to the ground state. 
In GaAs quantum dots Ti has a strong magnetic field dependence and can be very long, 
ranging from 1 ms around 5 T to more than 1 s at 1 T [123]. This dependence originates 
entirely from the spin-orbit interaction since, at such high values of the magnetic field, 
the hyperfine coupling plays no role for energy relaxation (due to the large mismatch 
between the nuclear and electron Zecman energies). At small magnetic fields the spin- 
orbit coupling becomes ineffective and, in fact, does not cause any relaxation at i? = 



(®) This value is a weighted average of the three nuclear species ®®Ga, ^^Ga, and '^^As, which 
have abundance 0.3, 0.2, and 0.5, respectively. For the three isotopes we have A — ^^/xsjU/Tjno, 
where m = (2.12, 2.56, 1.44) x /jjv, while ??Ga = 2.7 10^ and ??as = 4.5 10^ [122]. 
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[124]. Nevertheless, the hyperfine interaction contributes to the reduction of Ti to much 
smaller values (down to 10 — 100 ns, due to electron- nuclear flip-flops [40]). 

The transverse spin decay time T2 = = r~y describes the decay of the transverse 
polarization components and Py. The T2 time cannot be larger than 2Ti. This 
maximal value is obtained if only the spin-orbit coupling were present [124]. However, 
T2 is dominated by the hyperfine interaction and is much shorter than Ti. Due to the 
fluctuations of the Ovcrhauscr field in the nuclear bath's initial state, a transverse decay 
time of order 10 ns is obtained (see subsubsect. 4'2.1). In this case, it is clear that the 
much longer timescale tnuc = 10 — 100 /xs does not play a role for the transverse electron 
spin evolution. This decay time is usually denoted as T2 and referred to as 'ensemble- 
averaged' transverse spin decay time. We note that the decoherence process is generally 
non-exponential (see subsubsect. 4'2.1). 

If the initial nuclear state is prepared in an eigenstate of the Overhauser field in 
the B direction, an 'intrinsic' decay time T2 is obtained. A technique for narrowing 
the initial nuclear state was proposed in ref. [63] and is discussed in subsubsect. 4'2.3. 
The decay time T2 is determined in this case by the coupled dynamics of the electron 
spin and the nuclear bath. It is comparable to the tnuc time scale (estimates give T2 ~ 
1 — 100 lis) and therefore much longer than T2*. However, T2 is clearly very difficult 
to access experimentally. A quantity more easily measured is the spin echo decay time 
Techo- We refer to ref. [125] for a description of the spin echo technique, and to ref. [61] 
for its application to GaAs double dots. This method can be used to perfectly refocus an 
ensemble of spins in the idealized case where decoherence is only due to static fluctuations 
of the environment. However, in reality the initial polarization cannot be completely 
recovered due to the time evolution of the nuclear bath. A decay time Techo > 1 /US is 
reported in ref. [61] at 100 mT. 

3'4. Universal quantum gates. - Both single- and two-qubit gates have been demon- 
strated in GaAs quantum dots. The single gate was realized in ref. [126] by means of the 
well-known electron spin resonance (ESR), which we briefly describe here (for a more 
extended discussion see, e.g., ref. [125]). An oscillating magnetic field is applied in the 
transverse direction (perpendicular to B) at the resonant frequency uj = AEz/h. This 
ESR field can be seen as a sum of two contributions, rotating clockwise and counter- 
clockwise around B at the same frequency u. However, only the contribution precessing 
in resonance with the electron spin is of relevance. We denote this component by Bi, 
while the counter-propagating field is neglected in the following. 

Consider now the effect on the electron spin. Without the ESR signal, the spin simply 
precesses around B with angular frequency w. It is useful to introduce a reference frame 
rotating aroimd B in which the precessing spin appears static. We now apply Bi, which 
also appears static in the rotating frame. The effect is to induce a precession of the 
spin around Bi in the rotating frame. In particular, if the spin is initialized along B, a 
complete spin-flip is realized after a time -kH/ gfisBi. Typical fields in the experimental 
setup [126] are up to ^ 100 mT, which gives a ^ 15 ns switching time. 

On the other hand, the two-qubit SWAP operation was implemented with a much 



Quantum Computing with Electron Spins in Quantum Dots 



31 



faster gate duration ~ 0.5 ns [61]. The gate is realized in a similar way to the original 
proposal of ref. [7] (see subsect. 2'5), based on the control of the exchange coupling. In 
practice, the detuning e of the double dot is changed, since this modifies the splitting 
between the lower singlet branch and the triplet Tq, as described in subsubsect. 3 '2. 3 and 
illustrated in fig. 6. If a pulse from a large negative value eo to some value Eg around zero 
and back to so is applied, a finite energy splitting J(ec) between triplet Tq and singlet 
S exists for the duration r of the pulse. This causes the spin state to rotate to |J,t) 
if the pulse has length r = Trh/J{sc), which realizes the SWAP operation. The VSWAP 
operation is obtained if the time r is half of that required by the swap gate. 

3'4.1. Electrical manipulation of individual spins. While standard ESR is useful for 
single spin manipulation [126] and can in principle be applied to the individual dots of 
a large array (sec subsect. 2'5), it is much more convenient to perform coherent spin 
rotations through the electric gates at the individual dots. An example of such a tech- 
nique is the electric-dipole-induced spin resonance (EDSR), which is well known in two 
dimensions [127, 128, 129] and was also studied in lower dimensional-systems [130, 131]. 
EDSR in quantum dots was investigated theoretically in ref. [131], which discusses in 
detail the effect on the electron spin of an external ac electric field mediated by the 
spin-orbit interaction. In the following, we review the main results of this analysis. 

The single dot is described by the two-dimensional Hamiltonian 

(25) H=^ + U{r) + Hso + Hz - eEo • rsm{u;t), 

2m* 

where r = (x, y) is the electron's coordinate. The second term is the lateral confining 

potential of the dot and the third term is the spin-orbit coupling (Rashba and Dressel- 
haus) discussed in subsubsect. 3 '3.1. For the present section, it is convenient to define 
new axes Gx = (a^ + ay)/^/2, ey = — (a;^ — ay)/V2, and ez = a^ (instead of e, = aj, as 
in subsubsect. 3"3.1), where a^ are unit vectors along the cubic axes of the crystal. With 
this choice, Hso takes the particularly simple form 

(26) Hso = A ('-f' + ^ 

where A± are the spin-orbit lengths defined in eq. (22). The third term in eq. (25) is 

the usual Zceman coupling Hz = ff/^sB • S and the last term is the external electric 
perturbation. The electric field is assumed to be spatially uniform on the small region 
of the dot. 

The unperturbed states are the eigenstates ■0m(r)|±) of the dot Hamiltonian Hd = 
p^/2m* + f/(r). To calculate the effect of the oscillating electric field, one has to resort to 
third order perturbation theory since the final result has to be proportional to the spin- 
orbit coupling as well as to the electric field and to the Zeeman splitting /^Ez- No spin- 
electric coupling can be obtained at i? = 0, a property related to the invariancc of Hso 
upon time-reversal. It is convenient to approach the problem by making use of the unitary 
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SchriefFer- Wolff transformation He ^ described in tlie Appendix. If the confining 
potential is harmonic, i.e., U{r) = ^m*a;gr^, the final result is obtained explicitly as 
S = i£.-rT- z;#^(B X C) • where ^ = (t//A_, a;/A+, 0), C = (py/X-,p^/\+,0), and 
the spin Hamiltonian for the ground state reads 

(27) H^s = gfisB ■ S + gusho ■ S sm{Lot), 
where 

(28) ho = 2B X n = -^B X . 

Equation (27) clearly reveals the possibility to perform ESR-type spin manipulation, since 
the electric field induces an effective oscillating magnetic field ho sm{u>t). From the above 
expression, we estimate ho ~ 2 mT, using B = 1 T, A± = 10 /zm, and Eq — 100 V/cm. 
This value is in agreement with a recent experiment [132], in which a spin- flip time 
around 100 ns has been found. Corrections to the linear spin-orbit coupling and to the 
harmonic approximation of the confining potential are also considered in rcf. [131] and 
are responsible for an additional contribution to the EDSR signal which is not discussed 
here. This additional term only exists in combination with orbital effects of the magnetic 
field and is absent for an in-plane field B (as realized in ref. [132]). 

Finally, we note that EDSR is not the only method for spin manipulation via oscillat- 
ing electric fields. Spin-electric coupling can be also realized by oscillating the position 
of the dot in the presence of a static but inhomogeneous magnetic field, e.g., provided 
by the stray field of a nearby micromagnct [133]. This proposal was recently realized in 
ref. [134]. Gate- induced coherent single spin-rotations were also reported in [135] with a 
setup very similar to that of the EDSR experiment of ref. [132]. In that case however, the 
magnetic field is applied in the [110] direction, perpendicular to the electric field modu- 
lation, and eq. (28) gives ho = 0. Therefore, the spin-electric coupling in ref. [135] was 
attributed to the inhomogeneous Overhauser field produced by the hyperfine interaction. 

3'5. Readout of electron spin states. - Several methods are available for reading out 
the spin state of single and double quantum dots and all of them rely on the mechanism 
of spin-to-charge conversion. While the electron's magnetic moment is too small to be 
directly detected, the charge configuration of the single or double dot system can be 
measured accurately. This is usually accomplished by means of one or more quantum 
point contacts adjacent to the dots (these are narrow constrictions of the 2DEG through 
which current can flow), as shown in the right sample in fig. 4. The conductance of a 
point contact is quantized and is, at the transition between two plateaus, highly sensitive 
to the electrostatic environment, in particular to the charge distribution in the quantum 
dots (see the upper right panel in fig. 7). For example, a general strategy for spin 
readout consists in tuning the system to a configuration in which tunneling between 
different charge states is allowed or suppressed, depending on the particular spin state. 
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Fig. 7. - Spin readout in a single dot. In the upper left panels, two possible configurations 
of the dot are shown, which allow to distinguish between the spin states of the ground and 
excited levels. In (a) only an electron in the excited state can tunnel out, which is useful for 
the readout of Zeeman-split states in a single electron dot. In (b) one of the two states (here, 
the excited state) happens to have a much larger tunneling rate. This configuration is useful 
for discriminating singlet and triplet states in two-electron dots. In the upper right panel, the 
conductance of a point contact is shown as a function of the voltage of an adjacent gate. The 
cross marks the most sensitive position for charge detection. In the bottom panels, an example of 
the point contact current signal for energy selective readout is shown: The actual measurement 
is performed in the tread interval (marked in both panels by the horizontal threshold bar), during 
which the dot is as in (a). In the bottom right panel, a signal above threshold is detected as long 
as the dot is empty. This corresponds to the spin-||) electron tunneling out of the dot, followed 
by tunneling in of an ||) electron. In the left panel, conversely, the dot is always occupied by 
an electron in the ||) state. The point contact signal before tread refers to the initialization of 
the dot. See, e.g., ref. [40] for greater detail. (The top right and the two bottom figures are 
reprinted with permission from ref. [40]. Copyright (2007) by the American Physical Society.) 



The occurrence of the tunneling process is monitored by the point contact signal, and 
the corresponding spin state is inferred. We discuss below some more specific examples. 

3'5.1. Single dot readout. The first method for single-shot spin detection we will 
address here is the so called energy-selective readout [136]. The chemical potentials of the 
single dot are aligned with the 2DEG reservoir such that the |t) ground state is blocked 
while a ||) electron can still tunnel out (see fig. 7). Whether the electron leaves the dot or 
not is ascertained through charge sensing. It was proposed in [137, 138] to demonstrate 
the ESR rotation of a single spin by making use of this method. Note, however, that 
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the Zeeman splitting must be much larger than the thermal broadening, which implies 
large ESR excitation frequencies in the microwave range. This has caused additional 
problems, e.g., photon-assisted tunneling out of the dot, that are difficult to overcome 
experimentally and prevented this type of experiment from succeeding. Instead, the ESR 
experiment was performed in double dots using another type of spin readout, which we 
discuss in the next section. We would also like to mention that another detection method 
exists for single dots in the case where there is a large difference between the tunneling 
rates of two states [139]. This method is useful to discriminate singlet and triplet states 
in a two-electron dot {i.e., the 5(0, 2) and 7^(0, 2) discussed in subsubsect. 3'2.1), since 
the excited triplet has a more extended orbital wave function and better contact with 
the 2DEG reservoir. 

3'5.2. Spin readout in double dots. The mechanisms used for readout of double-dot 
states are identical to those described for initialization in subsubsects. 3'2.2 and 3'2.3. 
For example, if a double dot is in the charge configuration (1,1), tunnehng to S{0,2) 
is only possible if the spin state is |ti) or |J,t). Again, the (1,1) and (0,2) states are 
easily distinguished by charge sensing. This method of spin readout is used in the ESR 
experiment of ref. [126]: First, the double dot is initialized to a mixture of T+(l, 1) and 
T_(l, 1) via Pauli spin blockade (see left panel of fig. 5). The system is then brought 
to the Coulomb blockade regime by decreasing the detuning e. As a consequence, the 
chemical potential /is(0, 2) becomes higher than /Xa(l,l) and tunnehng is not allowed. 
The ESR signal can now be applied to one of the two electrons (^°), thereby rotating the 
initially parallel spin configuration (jtt) or to the antiparallel states jti) or |J,t). 
The double dot is finally brought again to the Pauli spin blockade regime and tunneling 
can now occur. Repeating the procedure many times, the probability that the single spin 
was rotated during the application of the ESR signal is determined. 

Another detection method can be used in the (1,1) charge configuration to distinguish 
jti) from Consider fig. 6, and suppose that at large negative detuning the system 
is in one of the two states jti) or ||t)- By adiabatically changing the detuning to large 
positive values, jti) evolves to 5(0,2) and |J,t) to To(l, 1) (for definiteness, we assume 
here and in the following ABn > 0) and the two different charge configurations can be 
distinguished. We can now describe the experiment realizing the VswAP operation [61] in 
greater detail: The system is first initialized in jtJ.) as described in subsubsect. 3'2.3. A 
pulse in e is then applied, which introduces a large singlet-triplet splitting, as discussed in 
subsect. 3'4. The swap or ^SWAP operations are realized for appropriate pulse lengths, 
but for an arbitrary pulse the double dot is brought, at large negative detuning, into a 
superposition of |t |) and ||t)- Finally, the detection method described above is applied 
and, repeating this scheme many times, the probabilities of the two spin states are 
measured. Singlet-triplet spin echo experiments can be also performed with a similar 
procedure [61]. 



(^'') Because of the hypcrfiric shift ABn of the magnetic field in the two dots, only one of the 
two electrons is usually in resonance with the ESR signal. 
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4. — RelcLxation and spin decoherence in GaAs quantum dots 

Electron spins in GaAs quantum dots arc inevitably coupled to the surrounding envi- 
ronment. This couphng results in decoherence, which is the process leading to the loss of 
information stored in a qubit. While an introduction to the mechanisms behind decoher- 
ence of electron spins in GaAs quantum dots has already been presented in subsect. 3'3, 
we will review this topic here in much greater detail. For the reader's convenience, we 
will initially repeat some of the basic concepts with additional remarks required for the 
subsequent treatment. 

As already mentioned in subsubscct. 3'3.3, two time scales describing the decoherence 
process of a single spin can be distinguished. The spin-flip time (also called longitudinal 
spin decay time) Ti describes the time scale for random spin-flips |t) whereas 
the transverse spin decay time T2 describes the decay of superpositions of spin-up and 
spin-down states a |t) +/3 | j). The first time scale Ti is important if the qubit is operated 
as a classical bit. For quantum computing, also the spin decoherence time T2 plays a 
major role and must thus be sufliciently long. 

The relaxation time Ti and the decoherence time T2 are not unrelated. Naively, 
one might expect T2 <C Ti, but as shown in [124] and as we will sec below, this is not 
necessarily the case for arbitrarily large Zeeman splittings. lu general, the electron spin 
S couples both to the external magnetic field B, and to the fluctuating internal field h(t) 
with the time-averaged value (h(t)) = (actual sources of this volatile internal field will 
be discussed later on). We also assume {hi{t)hj(t')) oc 5ij, where hi{t) with i = x,y,z are 
the components of the vector h{t). The Hamiltonian for the single-electron spin reads 



where g is the 5-factor (g = -0.44 in bulk GaAs) and hb = 9.27 x 10~^^ J/T is the 
magnetic moment of the single electron spin. The relaxation time can be expressed in 
the weak coupling limit as [140] 



where wz = guBB/h is the Zeeman frequency. On the other hand, the expression for 
the T2 time is 



Notably, the relaxation contribution {2Ti)~^ has been separated from the dephasing part 
incorporated in the integral. It was proven that for the spin-orbit interaction (discussed 

in the next section), up to linear order in momentum, h(t) -6 = and the effective 
magnetic field can only have fluctuations transverse to the applied B-field. As a result. 



(29) 



H = gnBS-B + S-h{t), 



(30) 




(31) 
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the integral yields zero and the upper bound on is realized, i.e., = 2Ti, independent 
of the origin of the fluctuations [124]. 

In contrast to the single-electron case, experiments performed on an ensemble of 
systems with different environments arc subject to additional dccohcrcncc. It is therefore 
required to introduce an ensemble-averaged transverse spin decay time T2, which is 
typically much shorter than the transverse single-spin decay time T2. One finds various 
other symbols for T2 in the literature, such as Tc (the correlation time) and Tm (the 
magnetization envelope decay time) , used to emphasize the non-exponential character of 
the decay. See subsubsect. 3'3.3 for a summary of theoretical estimates and experimental 
results for these various decoherence time scales. 

For quantum computers to work on a large scale, it is crucial to understand the 
microscopic mechanisms underlying dissipation and decoherence, and to devise efl'ective 
methods to reduce their impact on the electron's spin dynamics. This would allow to 
achieve longer coherence times and, in turn, reduce qubit errors. Two main sources 
of decoherence in GaAs arc to be identified, namely the spin-orbit and the hyperfine 
interaction, (i) The spin-orbit interaction couples the electron's spin to its orbital degrees 
of freedom. The orbital motion is influenced by lattice phonons, which provide a large 
dissipativc bosonic reservoir. In this way, an effective coupling between the electron's 
spin and the phonon bath is established leading to energy dissipation and decoherence. 
(ii) The Fermi contact hyperfine interaction couples the electron's spin directly to the 
surrounding bath of fluctuating nuclear spins. In a typical GaAs quantum dot the electron 
wave function overlaps with wave functions of approximately 10'' nuclei. The electron 
spin dynamics is thus strongly affected by the nuclear spin bath. These two decoherence 
mechanisms will be discussed extensively in the remaining part of this section. 

4'1. Spin- orbit interaction. - The main features of the spin-orbit interaction have 
already been discussed in subsubsect. 3'3.1. Here, we would like to emphasize that this 
interaction originates, in fact, from the relativistic Dirac equation and provides a direct 
coupling between the spin S and the momentum p. In vacuum, the spin-orbit term 
derived from the Dirac equation turns out to be 

he 

(32) if^^ = ___s.(px vy), 

where c is the speed of light, m is the free electron mass, e the electron charge, S = <t/2 
with (T being the vector of Pauli matrices, p is the canonical momentum, and V is an 
electric potential. Accounting for the lack of spatial inversion symmetry in the bulk GaAs 
crystal and assuming the presence of an asymmetric confining potential originating in 
the GaAs/AlGaAs heterostructure, the two-dimensional spin-orbit Hamiltonian reduces 
to the sum of the Rashba [141] and Dresselhaus [120] contributions 



(33) 



Hso = Hr + Hn = a {p^cTy - pyC^) + /3 {PxCTx - Py<yy) , 
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Fig. 8. - An electron in its orbital ground state, as determined by the two-dimensional confining 
potential (thick red curve). The probability to find the electron becomes negligible at distances 
greater than dot's diameter A ~ 100 nm. The ground state is separated from the first excited 



where a is tunable by external gates while /3 is a material constant stemming from 
the bulk inversion asymmetry. The tunability of a allows one, in principle, to achieve 
a = ±/3. In this particular case, Hso = ol {px T Py) {<^x ± o-y) and {ax ± Uy) is conserved. 
As a result, the electron's spin decouples from its momentum degrees of freedom. 

4' 1.1. Relaxation through phonons. We consider in this section the following two- 
dimensional model Hamiltonian for an electron in the conduction band 



Equation (34) above only differs from eq. (25) by the last term. Hd = p'^/2m* + U{r) 
describes the electron in the presence of a lateral confining potential, which is assumed 
parabolic, i.e., U{r) = ■m,*uj'^r'^ /2, and determines the ground state orbital wave function 
Tpoir) = exp(-r2/2A2)/AV7r, where A^^ ^ ^-i^(rn*a;o)2 + (eJ5^/2c)2. The probability 
density IV^olr)!^ to find the electron outside a circle of radius A/2 is negligible, where a 
typical dot diameter is A ^ 100 nm (see Figure 8). The second term of eq. (34) is the 
Zeeman Hamiltonian Hz — gfJ'B^ ■ S and the third one is the spin-orbit interaction Hgoi 
discussed briefly in the previous section. As discussed for eq. (25), it is convenient to 
express _ff in a rotated coordinate system, such that Hso takes the simpler form cq. (26). 
Note that, for a typical GaAs quantum dot, the spin orbit lengths A± ~ 1 — 10 /im are 
much larger than the diameter of the dot and the orbital level spacing h^/m*\^ ^ 
1 meV « 10 K far exceeds typical experimental temperatures /csT and Zeeman energies 
g^sB. The last contribution in eq. (34) takes into account two different types of electron- 



state by an energy gap ~ jm*}?' . 



(34) 



H — Hd + Hz + Hso + Hei-ph {t) ■ 
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phonon interactions [142] 



where the time-dependence is due to the phonon operators. In the above equation, 6, 
creates an acoustic phonon with wave vector q = (qy , Qz) , branch index j, and dispersion 
u)qj. Furthermore, pc is the sample density and F{qz) = \(p{z)\^e'"^'^dz, where ip{z) 
is the electron wave function in the direction z of the confinement (the full wave function 
is (p{z)ip{r), where r = {x,y)). Note that F(g^) equals unity for \qz\ <C and vanishes 
for \qz\ » d~^, where d is the size of the quantum well along the z axis. The couplings 
/3q,j are determined by the piezoelectric electron-phonon interaction as follows 

2tt 

(36) Pqj = — Z?"""?^ {q.enj (q) + g^e,,,- (q)) , 

Q K 

where is the electro-mechanical tensor, k is the dielectric constant, and e^,j(q) is 
the phonon polarization unit vector for branch j. For GaAs, = hi4 if the indices 

fii/rj are a cyclic permutation of xyz, where hi4 « 0.16 C/ni^. Otherwise, P^'^^ = [143]. 
Finally, the deformation potential electron-phonon interaction gives 

(37) Sqj = ^S'^'' (9^e^,j(q) + gi,e^j(q)) , 

where E^" is the deformation tensor. For GaAs one has S'*'' = 'BoSfj,^,, where Sq « 7 eV 
[143]. Therefore, the above expression simply becomes Sqj = So^j,i, since only the 
longitudinal branch j = 1 gives a non- vanishing contribution. 

As in subsubsect. 3'4.1, it is convenient to approach the perturbative treatment of the 
Hamiltonian (34) by making use of the unitary Schrieffer- Wolff transformation described 
in the Appendix. The final result of this procedure is an effective spin Hamiltonian of 
the form 

(38) //eff = (^o(r)l e^He-^ ]^o(r)) = gf^sB ■ S + gfxsS ■ SB{t) + ..., 

where spin-independent terms are omitted. In the above equation, (5B(t) = 2B x ^{t) 
where eq. (A.6) immediately gives^n(i) = (V'o(r)| [{{1- P)L-^^), Hei^ph{t)] |^o(r)). The 
vector ^ and the superoperators P and are defined in the Appendix. Note that the 
form of Hes above is the same as eq. (29). However, as pointed out in the beginning 
of sect. 4, there can be only transverse fluctuations of the effective magnetic field, i.e., 
S'B{t) • B = 0, to first order in the spin-orbit interaction. This property holds not only 
for phonons, but is valid regardless of the nature of the charge fluctuations, as is seen 
from the general form of eq. (A. 5). 



zqii-r 
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41.2. Energy relaxation. If the scattering events are not correlated, i.e., if the phonons 
emitted and absorbed by the electron leave the dot in a time Tc which satisfies d/s < 
Tc £ A/s (with s being the sound velocity), then the expectation value (S) obeys the 
Bloch equation 



(39) 



(S)=5MsBx(S)-r(S)+T. 



In this formula, the decay tensor T and the inhomogeneous term T can be derived in 
the Born-Markov approximation for a generic S'B{t) which fulfills (^B(t)) = 0, and are 
expressed in terms of the spectral function 



(40) 



Jij{uj) 



2 ? 

9 l^ B 

2h^ 



{5B,{{))6Bj{t))e-^Ht. 



In general, besides a term proportional to the spectral function J, the tensor T receives 
an additional contribution from elastic scattering of the electron spin. However, because 
of the transverse nature of the magnetic field fluctuations, this contribution vanishes 
identically. The final result for T\ and T2 reads 

1 2 

(41) — = — = R.e[Ja;a;(wz) + Jxx{-^z) + Jyyi'^z) + Jvy{-^z)\, 

-tl J-2 

where ujz = gfisB/h is the Zeeman frequency, as in eq. (30). The explicit expression for 
ReJxx(w) reads [124] 



(42) 



it/2 



(2A+m*w^) ^ 7rpcs| Jo 

^g-(wAsine)V2s' 



Fl i^cosl 

Si 



—2 



up— 2 



where N^^ = {e^^'^ — 1) ^ is the Bose distribution and sj is the velocity of sound for the 

branch j. In GaAs, the sj have values si » 4.7 x 10^ m/s and S2 



Furthermore, = ^o<^j,i with ^ 



V27r/i 



liK ^sin26'. 



_ *3 
~ 7 cV, = 3\/27r/(,i4K~^ sin 

Pse — 'dV^irhuK^^i'icos^ 6 — l)sin0 with hi4 



1 3.37 X 10^ m/s. 

2 6* cos 61, 02,0 = 
0.16 C/ni2 and 

K w 13. The result for Jyy{co) is obtained by substituting A+ A_ in the above 

expression. Here A± arc effective spin-orbit lengths given in [124]. For a magnetic field 
B in the z direction the simple result A± = X± is obtained, where X± are defined in 
eq. (22). Note that both Ja;^, and Jyy are multiplied by w^, thereby exhibiting a super- 
Ohinic behavior. 

By explicitly evaluating this rather cumbersome expression, a relaxation time Ti w 
(825 ± 275) fis is found at S = 8 T. This result is in very good agreement with the 
experimental value T^^^ = 800 /is at the same strength of the magnetic field [136] . The 
rather large uncertainty of the theoretical prediction is due to the measured value of the 
5- factor. 
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Fig. 9. - Left panel: The relaxation rate 1/Ti given by eq. 42 is plotted as a function of the 
magnetic field B. Contributions due to the piezoelectric effect with transverse and longitudinal 
phonons are plotted as dashed and dotted curves, respectively. The contribution due to the 
deformation potential mechanisms is shown by a dot-dashed curve. Right panel: The relaxation 
rate l/Ti is plotted for large magnetic fields, where the l/Ti oc behavior is suppressed. 



4'1.3. Magnetic field dependence of the relaxation rate. Another important result 
is obtained in the range of magnetic field below ~ 3 T, when the relaxation rate 1/Ti 
increases with the magnetic field as B^. This can be inferred from the presence in eq. (42) 
of the prefactor uj'^uj^ (note that, in the final expression for l/Ti, one has to set uj = loz). 
There are three contributions to this power-law behavior: (i) The fluctuating magnetic 
field is proportional to the external magnetic field, i.e., 5B^ oc B^. (ii) The phonon 
velocity is proportional to the phonon dispersion, i.e., Vph{w) cx uj'^ with u — gfisB. (iii) 
The spin-orbit Hamiltonian couples to the magnetic field through its momentum, i.e., 
Hso oa Pa oa B. The final result oc B^ was confirmed experimentally in ref. [123]. In the 
opposite limit, i.e. for magnetic fields larger than ~ 12 T, the power law is suppressed, 
since phonons are averaged to zero over the dot size = s/hfiBB <C A. These two 
regimes have been summarized in fig. 9. 

4' 1.4. Magnetic field angular dependence. For an arbitrary direction of the magnetic 
field B = B{sni9 cosip, sin 6* sin cos 0) it was found [124] that 

1 {a^ + I3'^){l + cos^ e) + 2a/3 sin^ sin 2(p 
¥i " /32Ti(0 = ^/2,a = O) ■ 

Noteworthy, there is an interference between the Rashba and Dresselhaus terms which 
leads to a diverging relaxation time Ti oo when a = P, 9 = tt/2, and ip = 37r/4. 
This result is valid to all orders in the spin-orbit interaction and is an effect of spin 
conservation, which occurs for the special condition a — f3 oi the spin-orbit couplings 
(see also the discussion of Hso, earlier in this section). 
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4'2. Hyperfine interaction. - As discussed in the previous section, the relaxation pro- 
cess of the spin polarization in GaAs quantum dots is dominated by the spin-orbit in- 
teraction, which couples the electron spin to the phonon bath. If no other effect were 
present, the upper bound for the decoherence time T2 — 2Ti would be satisfied [124]. 
Measurements of Ti reveal that ultra-long relaxation times ~ 1 s are achievable for mag- 
netic fields B ~ 1 T [123]. Hence, long decoherence times might be expected as well. 
Unfortunately, measured spin decoherence times arc considerably shorter and range from 
1 /US [61] to roughly 10 /xs [126, 144]. Thus, the spin decoherence in GaAs must be dom- 
inated by other effects. 

The major source of decoherence in GaAs was investigated theoretically in [55] and 
was attributed to the hyperfine interaction with the nuclear spins. In fact, this is the 
cause of decoherence in numerous candidate systems for quantum information process- 
ing applications such as quantum dots [145, 61, 126], Si:P donors [146], NV centers in 
diamond [15, 147], and molecular magnets [148, 149]. 

A possible way to limit this decoherence problem in GaAs might be to use holes instead 
of electrons [150, 151, 152], since the Fermi contact hyperfine interaction vanishes in this 
case. The detailed form of the hyperfine interaction for holes was recently studied in 
ref. [152] and receives contributions from the dipolc-dipole interaction and the coupling 
of the electron orbital angular momentum to the nuclear spins. It is indeed found to 
be smaller than for electrons (but still sizable) and of Ising type, differently from the 
electron's isotropic Hcnscnbcrg interaction. Another strategy against decoherence would 
be to employ materials such as C, Si, Ge and others, which do not host any nuclear 
magnetic moment. Finally, one can insist on GaAs, which is still the most common 
material for quantum information processing applications, and deal with the nuclear 
spins. In this section we follow this last approach and focus on the decoherence process 
caused by the hyperfine interaction. Based on its detailed understanding, we suggest 
possible schemes for its reduction or elimination. 

The following Hamiltonian describes the hyperfine interaction in a single quantum 
dot in the presence of a magnetic field applied in z direction B = Be^ and dipole-dipole 
interaction Hm between nuclear spins 

(44) H = h-S + bS, + eJ2l!+Hdd, 

i 

where h = A^I^ is the nuclear magnetic field also known as the Overhauscr field, Ai is 
the hyperfine coupling strength at site i, b = qhbB is the electron Zeeman splitting and 
e = qi^nB is the nuclear Zeeman splitting. The nuclear (Bohr) magneton is denoted by 
III {i^b) and nuclear (electron) gi-factor by gi (g). Typically, the electron orbital energy 
level spacing for lateral quantum dots containing one electron is much larger than the 
energy scale of the hyperfine interaction. Therefore, the electron remains in its orbital 
ground state and no orbital excitations due to the interactions with nuclei are possible. 

The hyperfine coupling is non- uniform since we have Ai oc )^/'(ri))^. For GaAs, the 
average coupling strength weighted by the natural abundance of each isotope (cf. subsub- 
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sect. 3'3.2) is ^ « 90 /ueV. The typical energy associated with the hyperfine interaction is 
then A/N k. 10® s~^, where N k, 10^ is the typical number of nuclei in contact with the 
electron spin. This scale has to be compared with the dipole-dipole interaction energy 
{{SHdd)^Y^^ ~ 10^ s^^, which is much smaller and will be neglected in the following. 

The nuclear Zeeman term can be formally eliminated from the Hamiltonian by trans- 
forming to a rotating reference frame [153]. One can also separate the longitudinal part 
Ho from the transverse flip-flop part V of the hyperfine Hamiltonian, which leads to the 
following expression 

(45) H' = Ho + V=ib' + K)S, + \{h+S- + h-S+), 

where h± = h^i: ihy and S± = i iSy are ladder operators for the nuclear field and 
the electron spin, respectively, and b' = b— e. From this representation, it becomes clear 
that the longitudinal component V describes fiip-fiop processes: If the electron spin and 
the spin of a neighboring nucleus are opposite, their direction can be simultaneously 
swapped. This flip-flop mechanism allows for electron mediated diffusion of the nuclear 
spins and determines fluctuations of the nuclear spins polarization on a timescale of order 
100 iJ,s. 

4'2.1. Ensemble averaged decoherence time. A central question to address is how the 
initial state of the nuclei affects the evolution of the electron spin. Assume that for times 
t <0 the electron spin and the nuclear bath are decoupled and described by the density 
operators ps(0) and p/(0) respectively. At t = 0, they are brought into contact over a 
switching time scale that is sufficiently short, e.g., the time taken to inject an electron 
into a quantum dot. The initial state of the entire system p(t = 0) is continuous at i = 
and, hence, p(0~) = p(0+) = ps (0)0/9/(0). Following ref. [153], we assume for simplicity 
nuclei with spin-j and consider three types of initial spin conflgurations, namely 

N 

(46) p\'\0) = \^j){^j\ with |V>,)=(g)(v^|tfc)+e**'=yr^|U)), 

(47) pf\o) = y: (^)fp (1 - /t)"-"^ livt) m , 

(48) pf \0) = \n){n\ with h,\n} =J2A,I^\n) = ^ \n) , 

k ^ 

where |tfc) and |J.fc) are the spin-up and spin-down eigenstates of the fc-th nucleus, 

determines the nuclear polarization p = 2/| — 1, is an arbitrary site-dependent phase, 
N is the total number of nuclei, and Nf is the number of the nuclei in the spin-up state. 
\N^) thus denotes any product state of the form |tiTT • • • ) with nuclear spins up and 
N — spins down, and A = J2k^k- 

We note that p')^\o) and pf\o) are both pure states but in the first one 1-0/) is 
chosen to render the z component of nuclear spin translationally invariant: \ \'ipi) = 
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(2/| — 1) /2 — p/2, while in the second one \n) is chosen to be an eigenstate of with 

(2) 

eigenvalue pA/2. On the other hand, p} (0) is a mixed state, which corresponds to 
an ensemble of product states where the N spins in each product state are selected from 
a bath of polarization p. 

Now we evaluate the nuclear spin dynamics under Hq, i.e., when the flip-flop term 
V can be neglected. This is justified, for instance, at large magnetic fields, because the 
energy gap between spin-up and spin-down states of the electron is much larger than 
the nuclear splitting between spin-up and spin-down states, thus making the flip-flop 
transition energetically forbidden. We also assume the simplified case of uniform coupling 
constants Ak = A/N. Since [Hq,Sz] = 0, {Sz)^ is constant. Instead, \Hq,S±\ ^ and 
the transverse components, {S±) ^ = {Sx)^ ± i {Sy)^, have a nontrivial time dependence, 
which can be evaluated by tracing out the electron and nuclear degrees of freedom (<S'±)j = 
Tr [e'^°*S'±e-*^°V(0)] • The final results read 

(49) {S^f^ = Y: (1 - f^f'""' e±<^'+^^(^^)/^^)*, 

(50) (5±)f = (5±)o e±<'''+f^/2)t^ 

where M{Nf) = 2Nf — N is the nuclear magnetization on a dot with N-^ nuclear spins 
up and we generally set h = 1 from now on. Note the similarity between randomly 
correlated pure states and mixed states, which yield the same final result eq. (49). It 
is also seen that no decay in eq. (50) is obtained for the eigenstate |n), but a finite 
transverse relaxation time follows from eq. (49), due to the average over the binomial 
distribution. One can further evaluate the expression in eq. (49), by direct application 
of the central limit theorem, and obtain 



This Gaussian decay occurs on a timescale ~ 5 ns, for a GaAs quantum dot with 
p2 < 1 and A'' = 10^. 

4'2.2. Electron Spin Resonance. To study the electron spin dynamics, the electron spin 

resonance (ESR) technique is very fruitful. In order to achieve the resonance condition, 
an alternating magnetic field Bac = B^^^x cos ujt is applied in the transverse direction, 
in addition to the static out-of-plane magnetic field B = Be^- The Hamiltonian of this 
system is 

(52) Hesr = {hz + b)Sz + bi cos{ujt)Sx, 

where hi = g^sBac- The first term is _ffo of eq. (45), but we neglected the nuclear Zecman 
splitting e. We have again assumed that spin-flip processes are not important since the 




44 



Robert Andrzej Zak, Beat Rothlisberger, Stefano Chesi, and Daniel Loss 




200 400 600 800 200 400 600 800 
Burst time (ns) Burst time (ns) 



Fig. 10. - Rabi oscillations for four different driving fields Bac The gray circles represent the 
experimental data obtained in the transport measurement as described in the main text. The 
best fit of eqs. (54) and (55) is represented by solid and dashed lines, respectively. Two fitting 
parameters were used: the saturation value and the phase shift. The latter turns out to be close 
to the theoretical value 7r/4. See ref. [144] for more information. 



static magnetic field is large. As discussed in the previous section, the nuclear bath is 
in a superposition of eigenstates \n) of the hz operator, with corresponding eigenvalues 
/i". The distribution of the ft" eigenvalues is Gaussian, with mean Kq — Trlp/Zi^} and 
variance a ~ A/^/N, defined by = TT{pj{hz — ho)^}. As in the previous section, pi 
is the density matrix of the nuclear system. 

The decay of the driven Rabi oscillations is found in a rotating-wave approximation 
(valid for (6i/6)^ <C 1) and is given by [63] 

(^3, P,,)4,C+/5:eo.(|< + l)+„(jI,) 



where C — 1/2 — y^y2^^b^/8aj exp {bl/8<j^) erfc [b i/VSa^j is a time-independent con- 
stant. The above formula holds for t ^ max {l/a, 1/bi, 6i/2cr^) and ho + b = w. Inter- 
estingly, the decay is a slow power law oc and exhibits a universal phase shift of 
7r/4. Even if the Rabi period exceeds the timescale r ~ 15 ns for the transverse spin 
decay, this result implies that Rabi oscillations are visible, due to the fact that the power 
law becomes valid after a short time t for bi sa a. The universal phase shift originates 
from the off-resonant contributions, which possess a higher Rabi period and shift the 
average oscillation in phase. These predictions, i.e., both the non-exponential decay and 
the universal phase shift, have been recently confirmed experimentally [144]. The elec- 
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tron spin state was detected by a transport measurement in a double dot configuration. 
Electrons are transported through the double quantum dot via transitions from the state 
with one electron in each dot to the singlet state with two electrons in the right dot, as 
extensively described in subsubsccts. 3'2.2 and 3'5.2. The states with even spin parity 
(parallel spins) block transport, whereas spins with odd spin parity (antiparallel spins) 
allow for transport. Given that the system is initialized in an even spin-parity state, the 
resonant oscillating transverse magnetic field rotates at least one of two electron spins 
and lifts the blockade. The time evolution of the probability to measure an odd parity 
state after a time t has been calculated for the following two cases. When both electrons 
are in resonance with the alternating magnetic field, one finds 



(54) 

bi 



1 + cos 



' ' ^3/2 



8aH 

Secondly, when only one electron is on resonance, one obtains 

The first result eq. (54) is valid for times t > max (l/a, l/6i, 6i/2o-^) ~ 20 ns for a 
1.4 mT nuclear field and bi < 2a. Note that the 1/t term, which becomes important for 
bi > (J, oscillates with twice the Rabi frequency. This is the result of the simultaneous 
rotation of both spins. The term which is dominant for bi < a, oscillates with the 

Rabi frequency and stems from the rotation of one spin only. As might be expected from 
these considerations, P^^^ contains exclusively the latter term. The comparison of these 
theoretical results with the experimental data is presented in fig. 10. 

It should be emphasized that the power-law decay and the universal, i.e., independent 
of all parameters, phase shift w/i are obtained with the nuclear field being static during 
a time much longer than the Rabi period. Therefore, the good agreement between the 
experiment [144] and theory confirms the assumption of a static nuclear bath. 

4'2.3. Narrowing of the nuclear state. The evolution of the electron spin governed by 
eq. (52) depends on the value of the nuclear field since the effective Zeeman splitting 
is given by ft" + b, where hz\n) = /i" |n). This means that the resonance condition 
6 + /i" — w = for ESR also depends on the nuclear field. Therefore, a measurement of 
the electron spin state determines ft" and the related state of the nuclei. 

The eigenvalues of the nuclear field in equilibrium obey a Gaussian distribution, as 
discussed in the previous section. That is, the diagonal elements of the nuclear spin 
density matrix pj (ft", 0) = (n| pi \n) = (-\/27ra) exp — (ft" — fto)^ /2a^ are Gaussian 
with mean fto and variance a. After the electron spin is initialized in a state jf) at time 
t = 0, the system evolves under the Hamiltonian Hesr until a measurement of the spin 
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is performed &t t = tm- Of special interest is the probability to find the electron spin in 
the orthogonal spin state ||) and a given nuclear eigenstate |n) 



where bi = gfisBi as defined before, and (5„ = b + h^ — uj is the measure of deviation from 
the resonance condition. On the other hand, the probability to find the electron spin in 
the state |J,) with arbitrary configuration of the nuclei is easily obtained by integrating 
out the nuclear field Pi{t) = J dh'^pi {h'^,0) P^{t). The measurement with outcome |J.) 
performed on the system results in the collapse of the diagonal part of the nuclear spin 
density matrix into 

(57) PI {K, 0) ^ {K, tm) = PI (K, 0) 



according to the basic rules of quantum mechanics. If the measurement has a time res- 
olution smaller than 1/bi {i.e., giving time-averaged values), the probability of outcome 
IJ.) as a function of the nuclear field eigenvalue hf^ reads P" = b\/2{b\ + 45^). In turn, 
the nuclear spin density matrix is multiplied by a Lorentzian with width bx and mean 
h'^ = CO — b. As a result, for bi < a, the nuclear spin distribution becomes narrowed and 
prolongation of the electron spin coherence is achieved. Analogously, if the measured 
outcome is |t), the Gaussian nuclear spin distribution is modified as follows: 

(58) PI {K, 0) ^ p(^'T) [K, tm) = PI iK, 0) ^_p^ l^2y 

Hence, the probability to match the resonance condition b + h" — uj — is considerably 
reduced. It has already been proven (of. subsubscct. 4'2.2) that the nuclear spin evolution 
is slow enough to allow for multiple measurements of the electron spin (each of them 
performed after re-initializion to |t)) over a timescale on which the nuclear spins may be 
considered static. Repeating the initialization and measurement scheme M times under 
the assumption of the static nuclear field, we arrive at 

(59) PI {h:,0) (K) = ^pi (/ij,0) (Pf )"^ (1 - PJ»)^-«^ , 

where a| denotes the number of times the state ||) was obtained. From the experimental 
point of view it should be easiest to narrow the nuclear field distribution by performing 
measurements with bi <C a. Provided that the electron spin was projected to ||), the 
narrowing has been achieved. Otherwise, the additional initialization-measurement cycles 
should be repeated until the narrowing is observed. The driving frequency should be 
adjusted after each initialization-measurement cycle to match the resonance condition in 
order to systematically move towards the narrowed state. Such an adaptive scheme is 
described in refs. [154] and [63]. 
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Fig. 11. - Average nuclear polarization in GaAs (p = ^ '^ii^i), where i runs over the A'" nuclei) 
as a function of temperature at different values of the magnetic field. 



4'2.4. Exact solution with a fully polarized nuclear bath. The idea of fully polarizing 
the nuclear system has attracted considerable attention, since in this special configuration 
the hyperfine-induced decoherence problem can be overcome. Furthermore, the system 
of an electron spin confined to a quantum dot with a completely polarized nuclear bath is 
solvable exactly [155, 156]. Despite the fact that the complete polarization of the nuclear 
spins is not accessible with current experimental methods (to reduce decoherence by an 
order of magnitude, polarization of above 99% is required [153], whereas the current 
world record is about 60% [51]), some schemes to achieve full polarization of the nuclei 
have been recently proposed. Figure 11 illustrates the temperature range that has to be 
accessed to directly polarize the nuclei in the presence of the external magnetic field [157]. 
Cooling to sub-mK temperatures is an experimental challenge which is currently actively 
pursued. At such low temperatures refs. [158] and [159] predict a phase transition to an 
ordered nuclear state, even in the absence of the external field. Of crucial importance 
is in this case the long range nature of the electron-mediated RKKY interaction, which 
is determined by nonanalytic corrections in the momentum dependence of the electronic 
spin-susceptibility [159, 160]. 

Consider a single electron confined to a quantum dot in its ground state. As discussed 
earlier, coupling of an external field B to the nuclear spins and direct dipole-dipole 
interactions among them are negligibly small. Hence, the Hamiltonian can be written 
as in eq. (44) with e = and Hdd = 0. Additionally, we assume that the system 
is initially in the state j^'o) = [JJ-TT •••) — \^) ^ ITTT ^-^-j the electron spin j-IJ.) 

is opposite to the fully polarized nuclear bath state [TTT ■ ' ')• Due to conservation of 
total angular momentum, the time evolution of the wave function is fully described by 
\^{t)) = a{t) l^-o) + J2k Pk{t) ItrtT • • ■ ife ■ ■ • ) with the time-dependent coefficients a[t) 
and I3k{t) obeying the normalization condition |a(i)]^ -I- |/?fe(i)|^ — 1- The second 
term in |^'(i)) is a coherent superposition of states with exactly one nuclear spin flipped 
and thus resembles a magnon excitation. 



We win study the correlator Co(<) = - (*o| 6S^{t)S^ = (l - \a{t)\'^) /2, which 
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describes the decay of the electron spin from its initial state j-U-). Here, 6Sz{t) = Sz{t)—Sz, 
where Sz{t) is the Heisenberg representation of Sz. Inserting |\t'(t)) into the Schrodinger 
equation, we obtain a set of coupled differential equations of the form 

(60) 2i^^=-Aa{t)/2 + Y,AkPk{t)-ba{t), 

k 

(61) 2i^^ = (A/2 - Ai) m + Aia{t) + 

with A = X^fe^fc ^^'■^ ^ = SfJ-sB. These equations can be solved in the usual way by 
performing a Laplace transform a{u) = dta{t)e~^*. One obtains 

(62) -W=*#^ + ^E^ ^'^'^'^ 



D{u) 2D{u) ^ iu- {A + 2b)/A + Ak/2' 



where 



(63) D{u)=iu+^^-\y-. rr—^TT. — r^- 

^ ' ^ ' ^ 4: 4^iu-(A + 26)/4 + Afc/2 

Defining iu = iiv + {A + 2b) /A, using the initial conditions a(0) = 1 and /3;c(0) = 0, and 
replacing the sum over k by an integral, i.e., 

(64) y 3^=2\A-2mNuj I dz\n(l-'-^^^\ 



where xo(-2) is the normalized transverse wave function defined through the envelope wave 
function |\l/(r, z)\^ = ex.\){—r'^/a?)x^{z)/{-Ka'^az), where a and are dot's dimensions in 
the lateral and transverse direction, respectively, one arrives at the following expression 
for a{t) (required to calculate the correlator Co(t)) 



(65) a{t) = — / djjj 



I 



Here, A' = A + 2b, and z is dimensionless in units of a^. The integration contour is a 
vertical line in the complex to plane such that all singularities lie to its left. There are two 
types of singularities: two branch points u) = and uiq = iAxo(0)/27rA'', and first order 
poles which lie on the imaginary axis {lu = itj). For 6 > there is one pole, while for 
b <Q there are two poles; for 6 = there is one pole at wi « iA/2 + iAj dzXQ{z)/4:nN. 
The contribution from the branch cut between u) = and u) = uq is 
(66) 

e-iA't/4 ri 2zoKe'-''^ 

<^{t) = 77 — / dn 5 , 

Jo |«/dz|x§(^)/xi(0)«-l| + «/7r7V-26Mxg(0)|' + (27r^o)'K2 
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where r' = tXo{0) with r = At/2TrN, and zq = zo{k) is defined through Xo{zo) = Xo(O)'^- 
N = noazO^ is a number of nuclear spins of density no and we have introduced the 
dimensionless variable k = lu/loq < I. 

The physical picture is the following. At time t = the system has some energy 
corresponding to the pole. It then starts oscillating back and forth, while each time 
visiting different frequencies within the branch cut. This corresponds to the flip-flop 
processes with the nuclei located at different sites. The contribution from the branch 
cut therefore describes the electron spin decoherence. At time r ~ 1 (where the decay 
mainly takes place) the decoherence is due to the interaction with the nuclei located 
at distances of order of the dot radius where the derivative of the coupling constant is 
maximal. For longer times, t ^ 1, the asymptotic behavior is determined either by the 
interaction with the nuclei located far from the dot or near the dot center depending on 
the Zeeman field value. 

First, consider the limit of large Zeeman field (5 ^ A) and large times (r S> 1). The 
main contribution to the integral comes from k —^ 1, i.e., by the interaction with the 
nuclei located near the dot center. The asymptotic behavior of a{t) is 



Here, (xo)" denotes the second derivative of Xo evaluated at ^; = 0. Remember that the 



the branch cut contribution a{t) (which is a decaying function of time), an oscillating 
term due to the pole contribution. Therefore, the leading term in the full correlator for 
r ^ 1 is a constant, given by the square of the modulus of the pole contribution. The 
next higher correction is the product of the pole contribution and a{t) as stated above. 

Secondly, for the magnetic field turned off {b = 0) the asymptotic behavior of the in- 
tegrand in eq. (66) for r 1 is determined by k <C 1. Taking for instance Xo(-^)/Xo(0) = 
, we find 



This result is non-universal. It depends on the form of the electron wave function at 
distances larger than the dot size, since the decoherence is due to the interaction with 
the nuclei located far from the dot. The disturbance of the nuclear spins propagates from 

the center of the dot outwards. 

To summarize this part, we have shown that the decaying part of the correlator Co{t) 
is strongly affected by the magnetic field strength. However, the characteristic time 
scale for the onset of the non-exponential decay is the same for all cases and is given by 
{A/N)~^. For GaAs quantum dots, this is in the range of microseconds. 

4'2.5. Time evolution with arbitrary nuclear polarization. In the previous section 
we discussed the exactly solvable case of an electron spin coupled to a fully polarized 



(67) 





contains, besides 



(68) 



a{T » 1) DC l/ln2 T. 
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nuclear bath. However, for an initial nuclear spin configuration which is not completely 
polarized, no exact solution exists and standard time-dependent perturbation theory fails 
[155, 156]. In this section, we review a systematic approach to the electron spin dynamics 
in the presence of the Fermi contact hypcrfine interaction eq. (45). This theory is valid 
in the limit of high magnetic fields and for arbitrary spin polarization and nuclear spin / 
(/ = 3/2 in GaAs). Furthermore, the initial state of the system is p(0) = ps(0) (S) Pi{0) 
and it is assumed that pi{0) = \n){n\, where \n) is an cigenstate of h^. In practice, this 
is not usually the case and an appropriate 'narrowing' procedure as the one described in 
subsubsect. 4'2.3 has to be applied. Throughout this section and in the next one we also 
assume isotropic hyperfine couplings (in d dimensions) 

(69) Ak = Ao exp 

where rk is the radial coordinate of the fc-th nuclear site. Note that J2k = A, and 
therefore Aq ~ A/N, where N is defined as the number of nuclei within the radius Iq and 
A ~ 90 /icV in GaAs. In the following, we choose the unit of energy such that Aq = 2. 

We follow ref. [153] and start from the generalized master equations describing the 
exact dynamics of the reduced electron spin polarization 

(70) {S,)t = N,{t) - i fdt'Jl.^it - t'){S,)t', 

Jo 

(71) {S+}t = icOn{S+}t - i I dt'i:++{t - t'){S+}t', 

Jo 

where a;„ = 6' + {n\hz\n) and S± = Sxi: iSy. It is also useful to introduce the reduced 
self-energy superoperator, which acts on a generic operator O as follows [153] 




(72) ±s{t)0 = -iTri {Le-'^-^*Ly[p/(0)O]} , 

where L = Lq + Ly is the fidl Liouvillian superoperator, defined by L{)0 = [iJo,0] 
and LyO = \y,0], the superoperator Q is defined by QO = O — p/(0)Tr/O, and Tr/ 
is the partial trace over the nuclear degrees of freedom. By using eq. (72) above, we 
can write the kernel of eq. (71) as Yj+^{t) = Trs[S'_S5(i)5'+], where Trg is the partial 
trace over the electron spin. We also define S^,y(t) = Tvs[p^T,s{t)Pu] where ii,v ^tji 
and = (1 ± (Jz)!'^- In terms of these quantities, the kernel of eq. (70) is given by 
Sz2(t) = 5]|i-(f) — Sfj^(f). Finally, the inhomogeneous term of eq. (70) is most simply 
defined in terms of its Laplace transform /(s) = dtf{t)e~^* , which gives Nz{s) = 

(2is)-MSTT(s) + 5^n(s)]- 

Solutions to eqs. (70) and (71) cannot be obtained exactly and, to proceed further, an 
expansion in powers of the flip-flop term V (see eq. (45)) is performed. The self-energy 
superoperator (72) can be expressed as a series in terms of Ly, where only even powers 
appear, i.e., Ss(t) = E^^(t)+E^^(t) + . ... At order 2(A:+1), the self-energy contribution 
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is suppressed at least by a factor A'^, where A = N/un becomes sniall at large values 
of the external magnetic field. Therefore, the Born approximation Ss(i) = T,g'{t) is 
justified in this limit. The time evolution can be determined from the Laplace transforms 
of eqs. (70) and (71), which involve the following lowest-order self-energies 

(73) Sf^^ (S) = -iNc+ [1+ (S - iuJn) +I-{S + icOn)] , 

(74) S{f (s) iNc^ [/_ (s - iujn) +I+{s + lUJn)] , 

(75) Si'i(s) = -iN [c_/+(s) + c+/_(s)] , 

where c± = /(/ + !)-(( m(m ± 1) )) and (( F(m) )) = Em=-/ P{m)F{m). Here, P(m) 
is the probability to find a nuclear spin with z-projection m and F{m,) is an arbitrary 
function (for example {{m)) = pi gives the polarization p of the initial state). The 
explicit expression of /± (s) reads 



1 At d x\\n 



(76) I±{s) = F = - / dx- 



I x\ 



where v = d/q — 1 and we performed the continuum limit using eq. (69). Two examples 
obtained inverting eq. (76) arc plotted in fig. 12 as functions of time. Note that t is in 
units of 2/Ao, which corresponds to a timescale of order N/A ~ 1 /xs. 

We describe next the case A <C 1, in which the solution has a simple form in terms 
of the small parameter S = N/u)^. The longitudinal component is given by {Sz)t = 
(5,)oo+«Tf^(i), where 



(77) {S,)^ = [1 -25Io{c+ + c_)](S,)o + 2pI5Io, 

(78) (jf "(t) = 25Re [{C+I-{t) + C-I+{t)) e"''^"*] . 

Here, /q = I±{t = 0), c± were defined before, and C± = c±{{Sz)q ± \)- The transverse 
solution is {S+)t = G°^''{t) + af^it), where 

(79) ar(t) = [1 -5h{c+ + c^)]{S+)oe'^-\ 

(80) cjf^t) = 5 [c+/_(t) + c.I+{t)] (5+)o. 

In both cases, the solution is expressed as a sum of a term which is constant in absolute 
value {{Sz)cc and C7°^'^(t) respectively), and a decaying part ((j^°'^(i) and a'^'^{t) respec- 
tively). The former is almost equal to the initial polarization, {Sz)o or (S'+)o, except 
a correction of order 6. Reintroducing dimensional units {i.e., multiplying by Aq/A), 
one obtains that the fraction of polarization which decays is 5 ^ A'^/cj^N. The time 
dependence of this decaying contribution is determined by the I±{t) functions. As illus- 
trated in fig. 12, the I±{t) can have a very different form, depending on the ratio d/q. 



52 Robert Andrzej Zak, Beat Rothlisberger, Stefano Chesi, and Daniel Loss 




t[2SMJ 



Fig. 12. - Left panel: Plot of Re[/+(t)//+(0)], as obtained from eq. (76). The solid line was 
calculated with d = 3 and m = 1, which is appropriate for a hydrogen s-type envelope wave 
function. The dashed curve refers to a two-dimensional Gaussian envelope wave function (d = 
m = 2). The asymptotic behaviors in the two cases are given in eqs. (82) and (81) respectively. 
Right panel: Geometrical factor f{d/q) defined by eq. (87), with d being the dimensionality and 
q is related to the shape of the ground-state wave function (see eq. (69)). 



In particular 

(81) /±(t»i)«^ for^<2, 

82 /± t > 1 oc -^-^ for - > 2, 

q 

where the first expression is determined by the nuclei close to the origin, and the second 
one from distant nuclei with a small coupling Ak- These two asymptotic forms are 
relevant in different physical situations, e.g., the first one for s-type hydrogenic functions 
of Si impurities {d — 3 and q = 1) and the second for a parabolic quantum dot in two 
dimensions {q = d = 2). Note that in eq. (81) the decay has an oscillatory character, 
diff'erently from eq. (82). 

Within the second-order Born approximation scheme discussed here, the time evolu- 
tion at generic values of the expansion parameter A < 1 offers a rich scenario, including 
terms with exponential decay, non-exponential decay, and undamped oscillations. The 
same method also reproduces the exact solution of the previous section, at p = 1, and 
can be extended to higher order in Ly. We refer to ref. [153] for the detailed treat- 
ment of these topics. As a final comment we note that the electron spin dynamics, as 
described at this level of approximation, is non-Markovian. For large times the electron 
spin polarization does not display an exponential decay but is constant in magnitude. 
The interesting part of the time dependence is given by the small decaying functions 
cT'^^^{t) and a'^'^{t), which also have a non-exponential charater, see eqs. (81) and (82). 
This picture is reexamined in the next section, where a different result (i.e., Markovian 
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decay) is found at very large time scales (T2 ~ uj'^N/A^, see also fig. 13). 

4'2.6. Markovian dynamics. At large magnetic fields, it is in principle possible to 
carry on the perturbative treatment of the previous section in terms of V to an arbitrary 
accuracy. Higher order corrections to the self-energy are multiplied by a small factor 
A and are expected to be small. However, it will become clear in the following that 
these initially small perturbative corrections can grow with time and become significant. 
Therefore, we reexamine the problem of the previous section (i.e., the time evolution of 
the spin polarization if the initial state \n) of the nuclear bath is an eigenstate of hz) with 
special attention to the time evolution at large time scales. Wc focus on the transverse 
component, generally described by eq. (71). Full decay of this polarization component at 
large times, of order WnN/A^ » N/A, is found in ref. [161] and an asymptotic dependence 
oc 1/t^ is derived in Ref [162]. It is therefore necessary to improve the treatment yielding 
eqs. (79) and (80) to describe the full decay of the transverse polarization. We outline 
here a different method of solution, which in fact predicts a Markovian spin decay at 
large times [163]. 

As a first step, consider the following second order effective Hamiltonian derived from 
a Schrieffer- Wolff transformation of eq. (44) 



where iv = b + and (/^) is the raising (lowering) ladder operator for the nuclear 
spins at site k. Differently from eq. (44), a site-dependent nuclear g'-factor gi^ is taken 
into account, such that 7fe = gi^f^N /gi-i-B, and the term Hdd is omitted. The Schrieffer- 
Wolff transformation was performed assuming that corrections to the diagonal part of 
Hes of order of ^ A'^/Nb can be neglected, but corrections of the same order to the non- 
diagonal part should bo retained. This is justified by the fact that the bath correlation 
time Tc ^ N/A <^ Nb/A^ is small compared to the time scale on which the diagonal 
corrections become relevant for b ^ A. 

In this effective Hamiltonian, the perturbing term is Ves = XS^, where X is given 
by X = l/2j2kjti^kAiIj^ Ij^ /u}. We can now proceed as in the previous section, and 
approximate the memory kernel I]++(t) of eq. (71) to lowest order in V^g . However, the 
result obtained in this case is more accurate: the second order Born approximation in 
Ves, applied to (83), retains contributions up to fourth order in the hyperfine couplings 
Ak- Differently from the previous section, the Markov approximation can be applied in 
this case, and leads to a non- vanishing decay. The decoherence rate 



is given in terms of the non-diagonal part dynamics X{t) = e '^^Xe*'^*. The expectation 
value is taken with respect to the initial 'narrowed' nuclear spin state and Aoj is a shift 



(83) 




(84) 
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of the precession frequency {i.e., = S.gH'^r.+Auj)*^ where Xt is a slowly varying 

envelope) which has to be determined self-consistently, see ref. [163]. Provided that the 
initial nuclear spin polarization is smooth on the scale of the electron wave function, the 
matrix elements of operators such as I^lj can be replaced by their average values and 
the correlator {X(t)X(0)) takes the form 

(85) {X{t)Xm = ^ ^4^fe-^(^-^0t, 

" k^i 

where c± arc as in the previous section and a>„ = (n|ti)|n) . We consider now an initially 
uniform unpolarized state with equal populations of all nuclear Zeeman levels. This 
means that ((m)) = and ((m^ )) = /(/ + l)/3, while a;„ = b. The final result takes 
the particularly simple form 

where 

(87) ^,'d\ _q flY'^-'Ti2d/q~l) 



qj d\3j (r(<i/g))3 ' 

which is valid for any dimensionality provided that d/q > 1/2. The geometrical factor 
dependence is shown in the right plot of fig. 12. In 3D it may describe a donor impurity, 
in 2D a lateral gated quantum dot, while in ID a nanotube or a nanowirc. Note that 
for q/d = 2 the decoherence rate diverges so that the decoherence time tends to zero, 
which signals the failure of the Makov approximation. Furthermore, 1 /T2 oc depends 
strongly on the nuclear spin so that systems with large magnetic moment like In (with 
/ = 9/2) exhibit faster decay than, for instance, GaAs (with / = 3/2). Finally, the 
validity condition for the Markov approximation, T2 > Tg ~ N/A is satisfied if A/b < 1, 
and coincides with the requirement for the Born approximation. 

In heteronuclear systems the decoherence time is given by the sum of decoherence rates 
weighted by the natural abundance of each isotope Ui squared, i.e., I/T2 = Y^^f^^i, if 
interspecies fiip-flops are neglected. The quadratic dependence on isotopic concentration 
is particularly striking. In spite of the fact that all isotopes in GaAs have the same 
nuclear spin and nominally the same hypcrfinc coupling constant, the decay is domi- 
nated by intraspecies flip-flops between As spins. This effect may both explain why 
only Ga spins have been seen to contribute to coherent effects in transport experiments 
through (In/Ga)As quantum dots [145] and why polarization seems to be transferred 
more efficiently from electron to As - rather than Ga - in GaAs quantum dots [164]. 

Finally, taking into account the discussion of the present section, we summarize in 
fig. 13 the functional dependence of the transverse spin polarization. There, four regimes 
are schematically depicted : (i) \i t < Tc ^ N/A, the power law decay discussed in 
the previous section apphes [155, 156, 153]. (ii) If Tc <C t <C t ~ bN/A^, a quadratic 
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Initial quadratic decay: 

Yao, Liu, Sham, PRB '06 



Fig. 13. - Schematic representation of the decay of the transverse polarization. The different 
regimes are discussed in subsubsects. 4'2.5 and 4'2.6. 



correction oc -[t/rf is obtained in ref. [161]. (iii) If r < t < [b/Af{N/A), the 

exponential decay oc e~*/"^^ discussed in this section is vahd [163]. (iv) If i ^ the 
long-time power law decay is oc [162]. 

Appendix 

The Schrieffer-Wolff transformation 

The Schrieffer-Wolff transformation is a very useful method to treat the coupling of 
external orbital perturbations to the spin states in the presence of spin-orbit interaction. 
It is applied to the theoretical treatment of EDSR in quantum dots (see Section 3'4.1) 
and spin relaxation due to phonons and other charge fluctuations (see Section 4). We 
give here some additional details on the subject (see also ref. [131]). 

The Hamiltonian treated here is similar to eq. (25) and (34). More specifically, 
we consider H = Hd + Hz + Hso + where the unperturbed dot Hamiltonian is 

Hd = p^/2m* -I- U{v) as before (the eigenstates are |n, ±) = ■0„(r)|±), with energy e„). 
The only difference is that V{t) refers here to a general orbital perturbation, whereas in 
eq. (25) and eq. (34) V{t) specifically is the external ac electric field and the electron- 
phonon coupling, respectively. 

Consider now a unitary transformation of H of the form 

(A.l) H = e^He-^ o:iHd + Hz + V{t) + [S, V{t)] + ... 

where S satisfies [Hd+Hz, S] = Hso, such that the spin-orbit perturbation Hso does not 
appear in H. The transformed Hamiltonian H is accurate to first order in the spin-orbit 
coupling, which determines the last term [S", y(t)]. The solution for S can be formally 

written in terms of the Liouvillian superoperators Ld(z) defined by Ld(z)-^ = [Hd{z), 
where A is & generic electron operator. Note that the spin-orbit coupling in the form of 
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eq. (26) can be written as Hso = i[Hd, iLacr ■ where ^ = {y/X-,x/X+,0) and 

S satisfies 

(A.2) (id + Lz)S = iLa<T ■ i- 

To every order in Lz, eq. (A.2) is solved by 

oc 

(A.3) S = {l-P)Y,{-L-^'LzTi<T-t 

n=0 

where L^^ = -i e*(^''+*'')*dt, with 0+, and P is such that {ruT\PA\n' a') = 

{na\A\n' a') if e„ = en' and {nu\P A\n' a') = otherwise. Therefore, PA commutes with 
Hi and Z/^P = 0. Furthermore, and Lz commute (since and Hz do) and 
and P commute with them as wch, fohowing directly from their definitions. By making 
use of these properties, one can verify that (A.2) is satisfied by (A.3). 

We can now calculate the effective potential [5, F(f)]. The first non- vanishing con- 
tribution to the commutator is obtained to first order in Lz and can be written in the 
following form 

(A.4) [S, V{t)] = [g^isa • (B x (1 - P)L-/i), V{t)] + . . . . 

Finally, we derive an effective spin Hamiltonian by calculating the expectation value 
(V'o(r)l-ff IV'o(r)) with respect to the ground state orbital wave function of H^. Neglecting 
spin-independent terms we obtain 

(A.5) Heff = ^B-<T + 5/XB(Bxn(t))-eT, 

where the first term is the usual Zeeman coupling and an additional effective magnetic 
field is produced by V{t) through n(t), which is defined by 

(A.6) m = (Vo(r)|[(l - P)L-,^i,V{t)]\M^)). 
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